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FOREWORD 


This  volume  contains  concise  versions  of  the  papers  presented 
at  the  Conference  on  Environmental  Effects  on  Antenna  Performance 
which  did  not  appear  in  Vol.  I  of  the  Proceedings.  The  papers  pre¬ 
sented  here  were  reproduced  from  the  manuscripts  supplied  by  the 
authors.  Editing  on  the  submitted  material  was  relatively  light, 
although  many  of  the  manuscripts  required  at  least  some  revision. 
As  editor,  I  should  like  to  thank  the  members  of  the  advisory  board 
for  their  assistance  and,  in  particular,  I  am  grateful  to  the  con¬ 
ference  secretary,  Mrs.  Eileen  Brackett,  for  her  diligence  in 
assembling  the  papers  and  readying  the  material  for  the  printer. 

Volume  I  which  was  issued  earlier  is  available  from  the 
Clearinghouse  for  Federal  Scientific  and  Technical  Information, 

Sills  Building,  5285  Port  Royal  Road,  Springfield,  Virginia  22151. 
When  ordering  specify  that  the  accession  number  is  PB-183  715. 
Additional  copies  of  this  Vol.  II  are  also  available  from  the  CFSTI 
for  $3.00  and  the  accession  number  is  AD  692  879.. 


James  R.  Wait 
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rROLATE  SPHEROIDAL  AND  LINEAR  AN  JENNAS  IN  LOSSY  MEDIA* 

R.  J.  Lytle 

Lawrence  Radiation  Laboratory,  University  of  California 
Livermore,  California 


Abstract.  The  short  circuit  current  and  charge  distributions  and 
the  complex  effective  length  for  a  prolate  spheroidal  receiving 
antenna,  in  a  complex  medium,  have  been  determined  for  an  arbitrary 
angle  of  incidence.  In  addition,  the  power  decay  as  a  function  of 
§  has  been  evaluated  for  a  variety  of ' situations  for  a  prolate 
spheroidal  transmitting  antenna. 


Tb  boundary  value  problem  depicted  in  Figure  1  has  been  solved 
by  utilizing  the  method  of  separation  of  variables. 


For  an  incident  wave  of  the  form 

(1)  =  E  e^k<^(-cos  +  sin  9. "a  ) 

O  X  X  1  z 

(2)  -ft1  .  ja 

R  'V 
c 


it  can  be  shown  that 


(3)  =  V  +  S  8  -  - 


J  E°  2  8<h 

r1  •  b  ""o — o — i  mn 

*1(s  -r)(s  -b 


(M)» 


*Work  performed  under  the  auspices  of  the  U.  S.  Atomic  Energy 
Commission. 
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r  p  3r^(c,§)  s 

<[(§  -l)  — -  cos  0  cos  m0  +  ro§  R^x'(c,§)  sin  0  sin  m0j  A  + 

L  °S  ran  mn 

?  dR^(c,§)  /.v 

Pran[(s  -l)  — ■— —  cos  0  cos  m0  +  m§  R^r;(c,§)  sin  0  sin  ra0]  - 

P  dR^(c,§)  /l  \  -) 

a  [(§  '-!)  — -  sin  0  sin  m0  +  m§  F,  ^(c,§)  cos  0  cos  m0]V-, 

mn  os  mn  j 

where  3  and  or  are  the  expansion  coefficients  for  the  scattered 
field,  and  is  4  jn  S^5  MJ/V^Jc) . 

As  E-j.  =  0  at  §  §o_,  by  multiplying  (3)  by  S^(c,T))  and  inte¬ 

grating  tne  results  from  -1  to  +1  in  i|  and  from  0  to  2tt  in  0,  and 
evaluating  the  result  at  §  =  one  can  obtain  for  Eq  =  k 


«  V?it 


+  /.n  [>  r -i 


-°  -  -i.  -,r°  )] 

J?-”  lt  ^ 


_  ro-~iSR-.4-;(c4c) 

lK-x  — ~~  +rf?r*Rit  (c^o)] 


w 


By  utilizing  this  result,  one  can  show  that  the  current 
distribution  1(11)  for  the  short  circuited  receiving  antenna  is 

-jlffTE  — 5 - - — 7~1 

(5)  isc(>l)  -  -nr2  J  (S^-i)'i-n2)  • 


(j)n  sW(=,\)  S^(c,I0 


(>*),„  ,  ,  .  «2.,,  “  In  H°' 


V  =  >  [!o  Ri„  <  =  ’U  +  K-V 


5  8  ‘!(=,U 


so 

The  short  circuit  charge  distribution  Q  (T))  can  be  determined  via  use 
of  the  continuity  equation  and  equation  (5). 
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By  appropriately  coupling  the  above  results  with  previously 
obtained  results  (Lyt le  and  Schultz, 1968)  for  the  current  distri¬ 
bution  1**00,  charge  distribution  Qf' r  (T| ) ,  and  input  impedance  2C 
of  a  prolate  spheroidal  transraitfinig  antenna,  one  car,  obtain  such 
quantities  as  the  complex  effective  length  and  the  current  and 
charge  distributions  of  a  loaded  receiving  antenna. 


Another  quantity  of  interest  is  the  rate  of  decay  in  the  com¬ 
plex  medium  of  the  power  generated  by  a  transmitting  antenna.  This 
has  been  determined  by  evaluating  the  power  passing  through  a  series 
of  confocal  prolate  spheroids  surrounding  the  transmitting  antenna. 


Examples  of  results  obtained  for  some  of  the  above  quantities 
follow . 


R.  J.  Lytle  and  F.  V.  Schultz  (1968),  Linear  Antennas  in  Plasma 
Media,  Purdue  University,  School  of  Electrical  Engineering. 
■TR-EE  68-10. 


Radiated  Power  Dependence  Upon  Riase  Length  and  Loss 

("a"  is  the  semiminor  axis  of  a  confocal  spheroid 

surrounding  a  transmitting  antenna  with  a  semiminor 

axis  "a  ") 
o 


■in  1ST  TJT  VT 

fit 


!>»)«  ■  V  <T‘4H  i  Wil  «T«  ■ 


.rtf  uf  i.is  fil  in 


Dependence  Upon  rfcas 
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NUMERICAL  STUDIES  OF  THE  EFFECTS  OF  NON- PLANAR 
LOCAL  TERRAIN  AND  GROUND  SCREENS 

by 

R.  V.  Row  and  D.  M.  Cunnold* 

Sylvania  Electronic  Systems,  Waltham,  Mass. 

ABSTRACT 

Numerical  integration  of  the  vector  form  of  Green's  theo¬ 
rem  for  electromagnetic  fields  has  been  used  to  evaluate 
the  coupling  between  a  pair  of  vertical  dipoles  over  a 
two-dimensional  non-planar  electrically  non-homogeneous 
ground.  Antennas  with  arbitrary  free-space  radiation  pat¬ 
terns  can  be  accommodated.  The  results  of  this  method  are 
compared  with  previous  computation  for  planar  ground  and 
ground  screen  using  the  compensation  theorem  approach. 
Agreement  is  good  except  at  angles  below  a  few  degrees. 

New  results  on  the  effects  of  a  finite  ground  screen  sys¬ 
tem  in  proximity  to  a  sloping  beach  and  a  flat  sea  with 
varying  tide  levels  are  given  and  compared  to  some  HF 
band  measurements. 


In  an  earlier  paper  (Cunnold,  Row,  Arnold,  1968]  the  wri¬ 
ters  described  a  method  of  determining  approximately  the 
coupling  between  two  horizontal  electric  (dipole)  antennas 
over  a  piecewise  homogeneous  non-planar  two-dimensional 
ground. 

The  method  if  based  upon  numerical  integration  of  Green's 
vector  theorem  for  electromagnetic  fields 

E  =  •  8J  { -iw/z(nxH)9  +  (nxE)  x  +  (n*E)Tty]ds  (1) 
e-ikR 

where  $  =  —  is  the  scalar  free-space  Green's  function 


*  Dr.  Cunnold  is  a  consultant  to  Sylvania 
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and  the  required  field  components  nxH,  nxE  and  n*E  are  ap¬ 
proximated  at  each  point  of  the  surface  in  the  sense  of 
classical  optics,  but  using  appropriate  Fresnel  reflection 
coefficients  at  each  point  of  the  surface.  Agreement  with 
a  limited  set  of  experimental  measurements  was  seen  to  be 
excellent.  Liepa  { 1968]  has  applied  a  similar  technique 
to  determine  the  field  from  a  horizontal  line  current 
source  illuminating  a  hemi-cylindrical  boss  on  an  infinite 
conducting  plane,  a  problem  for  which  an  exact  analytic 
solution  is  available.  He  found  excellent  agreement  be¬ 
tween  the  exact  and  approximate  field  distributions. 


The  approach  used  by  Cunnold  et  al  [  1968)  has  been  ex¬ 
tended  to  vertically  oriented  dipole  antennas  over  two- 
dimensional  ground,  as  shown  in  Figure  1.  As  in  the  ear¬ 
lier  work,  the  ground  surface  is  approximated  by  a  series 
of  connected  planar  strips  of  infinite  extent  in  the  y 
direction.  The  computer  program  evaluates  the  normalized 
transfer  impedance 


R  ~  R 
rad  rad 


(2) 


where  R  ,  is  the  radiation  resistance  of  an  elementary 
dipole,  0QZ§  is  the  transfer  impedance  between  the  anten¬ 
nas  in  the  absence  of  the  ground,  and  z£  is  the  ground 
scattered  contribution  from  the  i-th  strip. 

For  a  short  electric  dipole  of  length  df 

Z  =  E  df 
t  x 

Equation  (1)  is  used  to  evaluate  Ex  so  that  after  a  sta¬ 
tionary  phase  integration  in  the  y'  coordinate  only  the  z' 
integration  need  be  done  numerically. 

inn 

The  field  components,  Hyi,  E^, ,  E_(  and  Ex  are  calcu¬ 
lated  on  the  basis  of  Norton's  [  1937]  far-field  expres¬ 
sions. 
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Comparison  with  Known  Results 

Provided  the  Norton  ground  wave  terms  are  included  in  the 
field  expressions  ami  the  lower  antenna  of  the  pair  is 
taken  as  the  receiver*  excellent  agreement  is  obtained  be** 
tween  the  strip  method  and  Norton* s  results  for  dipoles 
over  a  homogeneous  flat  ground.  Wait  [1967]  has  examined 
the  situation  of  an  elevated  electric  dipole  over  a  cir¬ 
cular  ground  screen  on  a  flat  earth.  Agreement  with  Wait’s 
results  is  excellent  except  at  angles  below  5°,  where  the 
non-statlonary  nature  of  Eq.  (1)  will  exhibit  errors  in 
the  final  result  through  errors  in  approximating  the  true 
fields  nxH,  nxE  and  n*E  to  a  greater  degree  than  the  sta¬ 
tionary  formulation  used  by  Wait. 

Effects  of  Proximity  to  the  Sea 

An  antenna  installed  close  to  the  sea  is  in  an  environment 
which  changes  with  tide  level.  Figure  2  shows  a  profile 
for  which  extensive  computations  have  been  made  at  fre¬ 
quencies  of  5,  10*  20  and  35  MHz  and  screen  lengths  of 
150*  500  and  750  feet.  The  receiving  antenna  is  1  foot 
above  the  screen*  representing  a  monopole.  Figure  3  shows 
the  effect  of  Increasing  the  screen  length  from  150  feet 
to  500  feet  to  750  feet  under  low  tide  conditions.  A  sig¬ 
nificant  improvement  in  gain  is  realized  in  going  from  a 
150- foot  to  500-foot  screen*  with  the  maximum  improvement 
occurring  at  a  progressively  lower  angle  as  the  frequency 
is  increased.  A  lesser  additional  improvement  in  gain  is 
realized  by  extending  the  screen  over  the  down  sloping 
section  of  beach  to  the  edge  of  the  sea.  Figure  4  shows 
the  effects  of  a  10-foot  increase  in  tide  level  with  a 
150- foot  screen.  The  maximum  change  (increase)  in  gain 
is  about  1  dB  and  shows  up  at  the  highest  frequency.  Ad¬ 
ditional  computational  results  not  presented  here  show 
that  changes  in  gain  caused  by  tide  level  are  somewhat 
smaller  for  the  greater  ground  screen  lengths. 

The  computations  shown  in  Figures  3  and  4  were  repeated 
for  a  receiving  antenna  with  a  cosz  6  pattern.  Apart 
from  a  general  reduction  in  power  gain  at  angles  above 
about  10°*  the  results  are  in  substantial  agreement  with 
those  shown  in  these  which,  of  course*  pertain  to  a  dipole. 
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Comparison  with  Observation  of  Tide  Level  Changes 

Extensive  flight  test  pattern  measurements  [ Simons  et  al 
1968]  have  been  made  on  an  HF  band  log-periodic  monopole 
element  installed  on  a  low  profile  beach  installation  at 
Keflavik  Naval  Station,  Iceland.  The  contour  lines  are 
more  or  less  parallel  over  the  260°  to  340°  .seaward  sec¬ 
tor  and  hence  constitute  a  reasonable  approximation  to  a 
two-dimensional  ground  over  this  sector.  The  most  in¬ 
teresting  of  these  data  for  our  purposes  (shown  in  Figure 
7)  are  the  30-MHz  measurements  taken  on  October  29,  1967, 
at  an  azimuth  of  260°  (which  extends  over  the  sea).  One 
set  of  these  data  pertain  to  low  tide  and  the  other  to 
high  tide  (water  8-feet  higher).  Figure  5  shows  the 
actual  profile  ( solid  curve)  along  the  260°  azimuth, 
where  the  datum  le\el  is  transferred  from  mean  low  water 
(low  tide)  to  tne  14 -foot  contour  for  purposes  of  enter¬ 
ing  the  baseline  elevation  data  into  the  computer  pro¬ 
gram.  The  lightly  dashed  straight  line  is  the  sloping 
beach  profile  approximation  which  simplifies  the  labor 
of  entering  profile  data  into  the  computer.  The  strip 
method  program  was  run  for  both  these  profiles;  the 
former  (actual  profile)  case  being  designated  strip-by¬ 
strip  format  in  Figure  6.  In  the  computations  the 
ground  screen  (No.  10  AWG  wire  with  12-inch  spacing)  ex¬ 
tended  from  -55  feet  to  +80  feet.  It  is  evident  from 
Figure  6  that  the  results  obtained  with  the  sloping 
beach  format  are  in  excellent  agreement  with  the  strip- 
by-strip  input,  especially  for  the  low  tide  case.  The 
additional  set  of  curves,  shown  with  light  dash  and 
double  dots  in  Figure  6,  refer  to  the  gain  of  a  pair  of 
vertical  dipoles  6  feet  and  14  feet  above  an  infinite 
flat  sea,  respectively.  These  latter  curves  show  that 
in  general  it  is  not  possible  to  predict  the  influence 
of  changes  in  sea  level  by  such  a  simple  model.  There 
is  a  striking  resemblance  between  the  detailed  shape  of 
the  predicted  curves  as  shown  in  Figure  6  and  the  mea¬ 
surements  including  a  crossover  at  about  14°  in  both 
measurements  and  predictions.  The  maximum  increase  in 
gain  in  going  from  low  tide  to  high  tide  is  2.1  dB  at 
an  angle  of  10°  in  the  predictions,  and  3.6  dB  at  8.5° 
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in  the  measurements.  It  should  be  noted  that  at  angles 
below  1.5°  the  high  tide  measurements  show  a  relative 
field  strength  that  is  about  1.4  dB  greater  than  at  low 
tide.  The  opinion  is  offered  that  these  two  relative 
field  strength  (gain)  curves  should  be  coincident  at  all 
angles  below  lr5°r  as  is  the  case  for  dipoles  at  low 
heights  over  an  infinite  sea  shown  in  the  dash-double* 
dot  curves  of  Figure  6.  The  dash-double-dot  curve  in 
Figure  7  represents  a  uniform  downward  shift  of  1.4  dB 
for  the  high  tide  experimental  data.  When  this  is  done 
the  maximum  "observed"  increase  in  gain  in  going  from 
low  tide  to  high  tide  is  2.4  dB  at  9°,  which  compares 
even  more  favorably  with  the  strip  method  predictions  in 
Figure  6.  Of  course  the  "observed"  crossover  point  is 
now  shifted  down  to  12.5°  compared  to  about  14°  in  the 
predictions. 

It  is  concluded  that  the  strip  method  does  very  well  in 
predicting  the  detailed  shape  of  the  radiation  pattern 
at  angles  above  1.5°  or  2°  in  this  one  instance  of  a 
low-profile  beach  installation.  The  fact  that  the  mea¬ 
surements  were  made  using  the  LP  element  on  the  beach 
whereas  the  predictions  are  for  a  short  vertical  dipole 
should  not  change  the  conclusions  substantially,  as  evi¬ 
denced  by  the  discussion  in  the  preceding  section. 
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Figure  1.  Coordinate  Systems  for  Computing  Transfer 
Impedance 
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Figure  2.  Geometry  of  Merrimac  Light  Coast  Guard  Station  Sloping 
Beach  Profile 
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Figure  5.  Elevation  Profile  along  260°  Boresight  Azimuth  for 
Keflavik  Naval  Station 
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THE  CURRENT  DISTRIBUTION  ON  A  FINITE  LENGTH 
DIPOLE  IN  THE  PRESENCE  OF  GROUND  SCREENS 

Virgil  R.  Arens,  U.  R.  Embry,  and  D,  L,  Motte 

Sylvania  Electronic  Systems,  Western  Division 
Post  Office  Bex  205 
Mountain  View,  California  94040 

ABSTRACT 

Some  of  the  mathematical  and  numerical  techniques 
using  digital  computers  which  will  lead  to  determining  the  current 
distribution  on  a  finite  length  dipole  in  the  presence  of  radial 
ground  screen  are  presented.  Computed  results  for  the  input 
impedance  for  non-resonant  vertical  dipoles  near  tne  interface 
for  various  ground  conductivities  and  dielectric  constants  are 
given.  The  input  impedance  for  a  non-resonant  dipole  in  the  pres  - 
ence  of  a  planar  array  of  radial  wires  of  various  lengths  is  shown. 

L  INTRODUCTION 

On  a  number  of  occasions  investigators  (Wait,  1967) 
have  had  to  make  assumptions  with  respect  to  the  current  dis¬ 
tribution  on  finite  length  dipoles  above  and  near  variously  shaped 
ground  screens.  It  is  our  goal  to  develop  a  set  of  techniques 
that  will  allow  one  to  calculate,  using  digital  computers,  the 
current  distribution  under  these  conditions. 

The  methods  under  investigation  are  the  method  of 
moments  (Harrington,  1968),  Hallen's  integral  equation  (Hallen, 
1938;  Mei,  1965),  and  the  integral  equation  approach  (Tanner 
and  Andreas en,  1967).  All  three  of  these  approaches  make  use 
of  the  integral  of  free  space  Green's  function. 

.  2.  Green's  Function 

In  our  study  of  dipoles,  we  believe  that  we  have  devel¬ 
oped  some  useful  methods  in  evaluating  the  integral  of  the  Green's 
function. 

In  the  first  method  one  can  simply  perform  a  numerical 
integration  of  the  function.  Perhaps  one  should  not  use  the  word 
"simply",  though,  because  the  integration  for  the  diagonal  terms 
of  the  matrix  pose  some  extremely  difficult  numerical  problems. 
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For  these  terms  the  observation  point  lies  on  the  surface  of  the 
integration  segment  so  that  the  integrand  comes  within  5  of  a 
pole,  where  6  is  the  radius  of  the  wire.  To  get  an  idea  of  the 
difficulty  this  introduces,  the  approximated  numerical  value  for 
the  integral  has  not  become  stable  after  using  up  to  64  points  of 
Gaussian  quadrature  (corresponding  to  a  polynomial  of  degree  127  ) 
and  up  to  9000  points  ".sing  the  trapezoid  rule. 

In  the  second  method  one  rewrites  the  equations  in  cylin¬ 
drical  coordinates.  This  considerably  simplifies  the  work  because 
p  and  z  components  can  be  integrated  separately  after  a  local  rota¬ 
tion  and  translation  to  give  a  two  dimensional  problem  and  then 
dotted  with  a  direction  vector  at  the  observation  segment. 

The  third  method  makes  use  of  spline  functions  an11  has 
the  potential  of  allowing  us  to  develop  a  program  free  of  all  nu¬ 
merical  integrations.  This  is  accomplished  by  a  piece -wise 
approximation  of  e  -jkR  in  cubic  polynomials  over  one  cycle.  When 
the  resulting  polynomials  in  the  vc  riable  kR  are  placed  under  the 
integral  it  is  found  that  the  terms  ox  the  resulting  integrand  are  all 
integrable  in  closed  form.  The  particular  functions  one  uses  are 
called  cubic  splines,  an  excellent  account  of  which  is  found  in 
Handscomb  (1966). 

As  indicated  earlier,  when  the  observation  point  lies  on 
the  segment  of  the  cylinder  that  is  to  be  integrated,  a  singularity 
occurs  in  the  expression  for  the  electric  field.  In  point  of  fact  in 
the  derivation  of  the  equation  an  interchange  of  differentiation  and 
integration  was  made  which  is  not  necessarily  valid  for  these  seg¬ 
ments.  Fortunately,  YanBladel  (1961)  investigated  this  case  some 
years  ago  and  had  defined  a  principal  value  of  the  integral  for  the 
diagonal  terms. 


Evaluation  of  the  Integral  of  the  Green's  Function 


The  evaluation  of  the  integral  of  the  Green's  function  as 
was  indicated  earlier  can  at  times  be  evaluated  by  using  only 
numerical  techniques,  or  a  MacLaurin  series  expansion.  Some 
authors  have  in  the  past  tried  more  difficult  methods  that  involve, 
for  example,  Hankel  or  Bessel  functions.  Now  note  that: 
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With  this  simple  change,  the  integral  in  the  final  equation  can  now 
be  accurately  evaluated  by  as  few  as  two  or  three  points  of  Gaus  - 
sian  quadrature  for  all  values  of  R  and  accuracy  of  5  significant 
digits,  since  the  integrand  is  finite  for  R  equal  to  zero.  Here, 
then,  is  a  simple  and  accurate  numerical  technique  for  evaluating 
the  freespace  Green's  function,  which  we  have  not  seen  elsewhere 
in  the  literature.  It  illustrates  very  nicely  the  fact  that  it  is  not 
always  expedient  to  begin  a  numerical  evaluation  of  a  function  by 
expanding  it  in  a  power  series. 

4.  Antennas  Over  A  Finite  Conducting  Ground 

In  considering  antennas  over  ground,  it  is  easy  to  see 
that  the  free -space  Green's  function  is  no  longer  valid.  Here  one 
needs  to  match  a  set  of  boundary  conditions  at  the  interface  be  - 
tween  the  half  spaces.  A  Green's  function  which  must  meet  these 
conditions  is  called  a  modified  Green's  function.  Banos  (1966) 
gives  the  properly  modified  Green's  functions  that  are  valid  nearly 
everywhere  in  the  upper  half  space. 

The  total  Green's  function  is  the  sum  of  the  free  space 
Green's  function  for  the  image  of  the  source  segment  and  a 
correction  term  to  account  for  the  fi;;t  that  the  ground  has  finite 
conductivity.  Banos  deals  primarily  with  the  evaluation  of  this 
correction  term.  If  one  considers  a  dipole  just  above  the  ground, 
then  one  finds  that  none  of  the  expressions  given  by  Banos  are 
valid.  For  this  region  we  are  developing  an  adaptive  numerical 
integration  technique. 


5.  Ground  Screens 

Ground  screens  play  two  roles  in  high  frequency  antennas. 
The  first  role  is  used  in  association  with  monopoles  and  the 
ground  screen  then  is  commonly  called  a  ct  iterpoise.  Here  one 
is  interested  in  distributing  current  into  the  earth  in  a  least  loss 
manner  to  give  an  antenna  with  a  high  efficiency.  The  second 
role,  which  is  used  with  both  monopoles  and  dipoles,  is  the  use 
of  ground  screens  to  modify  or  change  the  far  field  pattern  of 
antennas. 

The  method  of  obtaining  the  current  distribution  in  the 
presence  of  ground  screens  can  be  obtained  in  at  least  two  ways. 

The  desired  method  is  dependent  on  the  number  of  wires  in  the 
ground  screen.  Consider  for  a  moment,  the  simplest  of  all  counter¬ 
poises,  a  limited  number  of  radial  wires  buried  in  the  earth  used  in 
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conjunction  with  a  monopole.  In  this  system  the  counterpoise  is 
considered  a  part  of  the  antenna.  The  impedance  matrix  which 
defines  the  antenna  includes  terms  representing  the  mutual  im¬ 
pedance  between  segments  of  the  antenna  and  the  buried  wires. 

When  the  number  of  wires  in  the  ground  is  large,  the 
size  of  the  matrix  may  become  too  large  to  consider  the  buried 
wires  as  part  of  the  radioing  system.  When  the  size  of  the 
ground  screen  is  sufficiently  large  one  can  take  the  expressions 
for  surface  impedance  given  by  Wait  (1956,  1958)  and  Larsen 
(1962),  and  determine  an  equivalent  dielectric  constant  and  earth 
conductivity  and  substitute  these  in  Bano's  equations  to  determine 
the  modified  Green's  function. 

If  the  ground  screen  is  not  large  or  the  edge  of  the  screen 
is  near  the  antenna  it  is  conceivable  that  one  could  determine  a 
different  modified  Green's  function.  However,  in  view  of  the  fact 
that  one  is  attempting  to  determine  the  mutual  impedance  between 
two  Hertzian  dipoles  due  to  the  presence  of  ground  screens. 

6.  Reducing  the  Number  of  Equations 

Use  of  the  symmetries  of  a  system  reduces  the  work  in¬ 
volved  in  solving  problems  connected  with  the  system.  In  simple 
structures,  the  symmetry  conditions  may  be  clear,  and  intuitive 
ideas  are  adequate  to  take  advantage  of  them.  In  a  more  compli¬ 
cated  system,  intuition  may  become  vague  and  a  systematic  method 
is  needed  to  examine  the  effects  of  the  various  symmetry  condi¬ 
tions.  We  have  developed  a  two  step  procedure  that  considers  all 
symmetry  conditions  when  setting  up  the  calculation  of  the  current 
at  specified  points  where  the  current  will  be  determined.  Symme¬ 
try  conditions  may  show  that  some  of  these  points  will  have  equal 
currents;  hence,  assign  only  one  variable  and  gather  the  coeffi¬ 
cients  to  form  a  single  term  in  each  equation.  (This  is  equivalent 
to  adding  columns  of  the  equation  matrix)  This  leads  to  an  over- 
deter mined  system,  with  more  equations  than  unknowns.  The 
second  step  is  to  determine  which  equations  are  redundant  and 
remove  them;  thus  obtaining  smaller  sets  of  equations  to  be  solved. 

As  an  example,  consider  a  monopole  with  radial  ground 
wires  spaced  ten  degrees  apart.  If  one  specifies  that  the  current 
will  be  found  at  five  points  on  each  wire,  the  result  is  a  set  of 
185  equations  in  the  185  unknown  currents.  If  the  symmetry  of  the 
ground  screen  is  invoked,  the  current  distribution  in  any  wire  must 
be  the  same  as  any  other,  so  five  unknowns  suffice  for  the  ground 
screen  current  values.  The  set  now  involves  ten  equations  in  ten 
unknowns. 
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7.  Results 


The  current  distributions  and  input  impedance  were  com¬ 
puted  for  a  center -fed -non -resonant  dipole.  The  dipole  length 
was  15  meters,  the  diameter  was  0.  006  meters  and  the  fre¬ 
quency  was  10  mHz.  For  compactness,  only  imput  impedance  plots 
are  shown*  For  most  cases,  the  current  distribution  obtained  was 
not  sinusoidal. 

Figure  1  shows  the  input  impedance  oi  the  dipole  in  free 
space,  and  with  the  tip  0,  01  meter  above  the  surface,  in  a  vertical 
position,  over  various  types  erf  ground.  The  value  fpr  the  perfect 
conductor  was  checked  by  computing  the  impedance  for  each  of  two 
collinear  antennas  fed  in  phase. 

Figure  2  shows  the  effect  on  the  input  impedance  when  a 
radial  array  of  wires  is  placed  near  one  end  of  the  dipole.  Values 
are  shown  for  various  lengths  of  4-wire  and  36-wire  arrays.  The 
values  for  9-wire  and  18-wire  arrays  show  similar  patterns.  As 
might  be  expected,  many  short  wires  have  more  effect  than  a  few 
short  wires.  As  the  wires  in  the  array  approach  a  length  of  half 
a  wave-length,  the  arrays  with  few  wires  cause  a  greater  change 
than  the  arrays  with  many  wires. 

Placing  a  vertical  dipole  over  Bmooth  homogenous  ground 
appears  to  increase  the  input  resistance.  Placing  the  dipole  over 
an  array  of  short  radial  wires  increases  the  reactance,  with  a 
small  change  in  the  resistance.  The  combined  effects  are  ex¬ 
pected  to  cause  the  input  impedance  to  approach  that  of  the  dipole 
over  a  perfect  reflecting  surface.  If  the  array  contains  many 
wires,  the  length  should  be  about  0. 1  5  to  0.  2  wave-length;  with 
few  wires,  the  length  should  be  greater,  to  about  0.  25  wave  length. 

9.  Conclusion 

The  digital  computer  programs  that  are  being  developed 
will,  when  complete,  allow  one  to  evaluate  the  properties  of  any 
dipole  antenna  in  any  orientation  between  horizontal  and  vertical 
over  a  ground  screen. 

The  approaches  are  being  developed  with  the  idea  in  mind 
that  other  antennas  such  as  the  log  periodic  or  rhombic  are  more 
important  than  dipoles;  the  work  on  dipoles  is  a  means  to  develop 
the  techniques. 

Besides  the  advantage  of  speed,  economy,  and  flexibility 
relative  to  experimental  methods,  the  methods  outlined  here  will 
provide  answers  with  a  precision  and  detail  unobtainable  by  flight 
testing  full  scale  high  frequency  antennas.  The  knowledge  of  the 
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current  distribution  and  input  impedance  provides  insight  to  the 
antenna's  operation  and  enhances  the  engineer's  design  tools. 
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MEASURED  PATTERNS  OF  HF  ANTENNAS  AND 
CORRELATION  WITH  SURROUNDING  TERRAIN 

D.  R.  McCoy,  R.  D.  Wengenroth 
General  Electric  Company,  HMES,  Syracuse,  N,  Y. 

and 

J.  J.  Simons,  RADC,  Griff iss  Air  Force  Base,  Rome,  N.  Y. 

A  large  number  of  fullscale  radiation  patterns  have  been  measured  on  several  HF 
antennas  with  various  siting  conditions.  From  the  data  it  can  be  demonstrated  conclusively 
that  siting  plays  a  major  role  in  the  radiation  characteristics,  as  expected.  A  study  of  the 
data  reveals  ways  of  modeling  the  earth  so  as  to  adeqvately  predict  the  influence  of  various 
parameters  on  antenna  coverage.  Thus,  these  data  provide  a  firm  reference  with  which 
various  theories  can  be  verified.  Only  two  examples  from  the  many  available  are  presented 
in  this  paper. 


1.  Introduction 

It  has  been  well  established  by  theory  that  antenna  surroundings  will  modify  the  far- field  radiation  patterns 
of  an  antenna  system -that  in  fact,  the  antenna  is  composed  of  the  actual  radiators  and  any  territory  many  miles 
away  which  may  be  illuminated  by  the  system.  In  the  high  frequency  range,  where  the  wavelengths  are  quite 
large  and  the  antenna  beam  widths  in  general  relatively  broad,  this  phenomena  becomes  extensive.  This  paper 
presents  some  results  of  a  measurement  program  on  such  a  system,  demonstrating  the  effect.  Only  a  brief 
sampling  of  the  data  is  presented  here.  A  more  complete  presentation  is  made  elsewhere/2,3)  it  also  discusses 
techniques  whereby  predictions  may  be  made  of  the  impact  of  surroundings  on  such  a  proposed  installation. 

2.  Discussion 

The  full-scale  patterns  were  measured  using  a  KC-135  aircraft,  instrumented  with  a  receiving  antenna  on  the 
tail  boom,  and  a  set  of  four  receivers  multiplexed  inside.  Four  transmitters  operating  at  different  frequencies 
were  time-multiplexed  into  the  antenna  on  the  ground,  allowing  for  four  patterns  to  be  measured  on  each  pass 
over  the  antenna.  The  aircraft  was  flown  at  a  constant  altitude  (nominally  30,000  ft.)  and  the  range  and  signal 
strength  were  recorded  in  the  air.  Flight  path  was  controlled  by  ground  radar  and  the  navigator  out  to  a 
maximum  tange  of  approximately  2$0  nautical  miles.  The  data  was  then  scaled,  calibration  data  applied  and  the 
elevation  angle  and  signal  strength  calculated  and  plotted,  yielding  relative  data  which  in  most  cases  was 
reproducible  to  within  1  dB. 

These  full-scale  patterns  of  HF  antennas  instatted  with  various  siting  conditions  have  been  made  jointly  by 
General  Electric  Company  and  the  RADC  test  center  The  installations  range  from  a  site  atop  a  mountain 
overlooking  the  ocean  in  Greenland  to  a  sea  shore  site  in  Iceland,  with  the  ground  screen  end  touching  the  ocean 
water  at  high  tide.  Most  of  the  results  discussed  here  are  of  two  identical  sets  of  log  periodic  elements,  one 
looking  north  and  the  other  south  over  land  in  upper  New  York  State.  Of  particular  interest  is  the  series  of 
patterns  obtained  at  boresight  (188°)  and  at  plus  and  minus  10°  (198°,  178°)  of  boresight  looking  south.  The 
antenna  is  on  a  hill  overlooking  the  Mohawk  valley.  From  the  elevation  of  1760  feet  above  sea  level  the  land 
drops  off  rapidly  to  7-900  feet  (a  slope  of  approximately  6°)  into  the  valley  10  to  20  miles  away.  In  the  198° 
direction  the  terrain  is  broken  up  with  several  hills  so  that  there  is  no  flat  land.  This  pattern  (not  shown)  does 
not  display  the  rapid  lobing  at  the  higher  angles.  The  profile  shown  in  figure  1  demonstrates  the  condition  on 
boresight  (188°).  The  radiation  pattern  for  this  direction  is  shown  in  figure  2.  The  identical  model  of  antenna 
looking  to  the  north  from  the  same  hill  exhibits  entirely  different  radiation  characteristics,  as  demonstrated  in 
figure  3.  In  this  direction  the  initial  slope  is  about  6°  but  the  elevation  levels  off  at  about  1400  to  1500  feet  in  a 
short  distance  and  stays  constant  to  the  Adirondack  Mountain  region  30-40  miles  away.  This  contouring  is 
shown  in  figure  4. 

A  consideration  of  the  specular  reflection  from  a  combination  of  the  slope  in  close  and  the  relatively  flat 
terrain  further  out  leads  to  a  rapidly  varying  pattern  at  low  angles  and  a  wide  minimum  at  slightly  higher 
elevation  angles.  It  therefore  appears  from  an  observation  of  measured  data  that  there  is  a  third  effect  at  the 
lower  angles,  in  which  the  flat  plane  ripple  is  suppressed.  This  is  probably  due  to  the  energy  being  blocked  by  the 
slight  roll  immediately  in  front  of  the  antenna.  There  were  minor  differences  in  the  lobing  structure,  particularly 
to  the  south,  when  snow  was  on  the  ground.  Several  values  of  ground  constants  were  considered  in  calculating 
the  above  mentioned  patterns.  Proper  choice  of  ground  constants  (K  =  30  and  O  -  .001)  seemed  to  give  a  good 
fit  to  the  measured  data,  with  the  reservations  of  suppression  of  rapid  lobing  at  the  lower  angles. 
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Figure  5  shows  the  measured  patterns  at  the  KcOavik  (iceiand)  site.  The  antenna  is  located  on  the  shore  line 
with  the  ground  plane  extending  to  the  water.  In  this  ease  there  was  a  noticeable  effect  of  the  tidal  variations, 
which  were  greater  than  13  feet.  A  calculation  of  the  expected  patterns  for  the  low  and  high  tide  conditions 
provides  an  excellent  prediction  of  the  effect  as  can  be  seen  from  figure  6.  This  same  antenna,  when  set  on  a 
2560  foot  hill  overlooking  the  snow,  rock,  hills  etc.  and  finally  the  ocean  17  miles  away,  at  Thule,  Greenland, 
displayed  violent  lobing,  even  though  the  earth  dropped  off  sharply  in  frc..t  of  the  antenna  to  an  elevation  of 
approximately  1500  feet  above  sea  level. 
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A  technique  of  modeling  the  land  in  front  of  an  antenna  with  a  scries  of  flat  planes  and  knife  edges  has  been 
developed/  0  This  program  estimates  the  worst  case  limit  for  the  particular  geometry  assumed.  A  result  of  using 
this  program  to  analyze  the  Thule  antenna  is  shown  in  figure  7.  In  this  case  the  very  rough  terrain  was  assumed 
to  be  approximated  by  two  knife-edge  points  and  intervening  flat  planes.  As  can  be  seen,  the  low  angle  coverage 
prediction  is  quite  good.  It  is  obvious  that  an  added  effect  occurs  in  the  2.5°  to  4.5°  region,  causing  a  wide  lobe 
of  high  amplitude  and  phasing  with  the  more  regular  pattern.  Perhaps  a  superposition  of  these  data  with  returns 
from  additional  specular  points  would  be  in  order  to  further  refine  this  result. 

In  the  case  of  the  data  presented  above  for  the  Stan  Hill  site  in  New  York  State,  there  are  no  significant 
sharp  obstructions  in  the  path  to  the  north  and  the  program  does  not  apply;  however,  it  has  been  used  with  fair 
success  on  the  south-looking  antenna/2)  The  program  utility  can  be  enhanced  by  more  accurate  modeling  and 
perhaps  by  replacing  the  knife  edge  diffraction  with  hemicylinders  and  allowing  for  extended  specular  reflection 
planes.  This  is  being  studied  at  the  present  time. 
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Figurt  7 


3.  Conclusions 


From  the  observed  data  it  is  completely  demonstrated  that  siting  plays  a  major  part  in  determining  the 
radiation  characteristics  of  an  antenna  working  in  at  the  HF  frequency  range.  It  is  also  evident  that,  by  carefully 
studying  the  terrain  condition  and  contour  for  several  miles  in  front  of  the  antenna,  a  good  prediction  of  the 
performance  can  be  made. 
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ABSTRACT 


This  paper  describes  the  performance  of  a  VHF  dipole  antenna 
located  at  the  tip  of  an  otherwise  pointed  reentry  nosecone.  Trans¬ 
mission  and  reflection  measurements  made  during  a  recent  flight  ex¬ 
periment  are  presented  and  compared  with  calculations  made  for  an 
appropriate  theoretical  model. 


INTRODUCTION 

The  performance  of  a  thin  (k0a  =  .092)  VHF  dipole-type  antenna 
located  at  the  tip  of  a  9°  nose  cone  during  an  almost  vertical  re¬ 
entry  at  a  relatively  constant  velocity  of  16.8  kft/sec  is  described. 
The  reentry  period  reported  covers  a  time  span  of  eight  seconds  with 
an  altitude  range  of  230  to  100  kft. 

The  antenna  as  shown  in  Fig.  1  is  an  asymmetrical  dipole  with 
the  cylindrical  quarter-wave  stub  acting  as  one  arm  and  the  payload 
skin  as  the  other  arm.  It  is  bal un-fed,  the  outer  conductor  of  the 
cable  being  connected  to  the  cone  and  the  center  conductor  to  the 
cylinder.  A  quarter-wave  choke  section  was  introduced  behind  the 
conical  arm,  hopefully  to  decouple  the  aft  portion  of  the  conical 
payload.  That  this  was  not  fully  achieved  is  evident  from  the 
free-space  radiation  pattern  of  rig.  2.  Note  the  "notch"  at  broad¬ 
side  which  is  primarily  attributable  to  the  excitation  of  the 
afterbody. 


*The  work  reported  in  this  abstract  was  supported  in  part  by  Contract 
F33615-67-C-1643  between  Air  Force  Avionics  Laboratory,  Air  Force 
Systems  Command,  Wright-Patterson  Air  Force  Base,  Ohio  and  The  Ohio 
State  University  Research  Foundation. 
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THE  FLIGHT  DATA 

The  altitude,  velocity  and  angular  coordinates  of  the  trajectory 
were  obtained  by  four  tracking  radars,  2750  to  5500  MHz.  The  X- 
band  telemetry  signals  were  recorded  by  three  receiving  stations. 

The  VHF  transmission  was  received  by  four  slightly  separated  CP 
antennas  and  provided  separate  VHF-AGC  records. 

The  telemetered  signals  of  an  on-board  gyro  and  accelerometer 
together  with  known  initial  reentry  conditions  and  the  precession 
frequency  provided  angle-of-attack  data  so  that  the  aspect  angle 
(the  angle  between  the  cone  axis  and  receiver  direction)  could  be 
computed.  The  ^receseion  rate  was  verified  by  the  cyclic  variation 
of  the  VHF  transnission  level  due  to  the  "notch"  in  the  pattern. 

During  the  flight  the  impedance  of  the  antenna  was  measured 
with  a  four-probe  refiectometer  and  telemetry  recorded  as  shewn 
in  Fig.  3.  Except  for  one  short  time  interval,  preliminary  to 
maximum  signal  attenuation,  the  modulus  of  the  reflection  coeffi¬ 
cient  was  relatively  constant.  In  addition  the  antenna  was  operated 
in  a  receiving  mode,  a  radiometer  being  used  to  record  the  effective 
noise  temperature  of  the  plasma  at  the  antenna  terminals.  This  data 
is  presented  in  Fig.  4  where  the  "average"  noise  temperature  is 
approximately  3500OK.  Outside  of  the  time  period  shown  the  antenna 
receiver  was  saturated  by  VHF  interference,  so  at  the  time  the  VHF 
transmission  was  undergoing  maximum  attenuation  the  plasma  noise  and 
interfering  signal  contribution  could  not  be  resolved. 

In  Fig.  5  the  VHF  transmitted  signal  which  was  recovered  from 
the  calibrated  AGC  records  of  three  different  receivers  is  pre¬ 
sented.  The  solid  curve  connects  the  average  received  values.  A 
maximum  signal  attenuation  of  17  dB  is  observed  at  193  kft. 


THE  ANALYTICAL  MODEL  OF  THE  ANTENNA 


We  approximate  the  flight  situation  with  a  gap-excited  finite 
dipole  in  an  infinite,  uniform  plasma  column.  The  far-zone  result 
is  similar  to  that  for  the  Hertzian  dipole  in  the  infinite  column, 
with  the  exception  of  the  current-plasma  interaction.  Using  the 
result  of  Ting*  et  al .  The  ratio  of  the  power  with  and  without 
plasma  radiated  in  the  direction  e  is  given  by 
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Fig.  >•  THE  antenna  and  payload. 


Fig.  2.  e- PLANE  RADIATION  pattern 
OF  antenna  (  f  *  231.4  MHz  ). 


Fig.  3.  ANTENNA  'MPEOANCE 
DURING  REENTRY. 


Fig.  4.  EFFECTIVE  PLASMA  NOISE  TEMPERATURE 
AT  ANTENNA  TERMINALS. 


where 


(2) 


1  +  i  v/oj 


2L  is  the  dipole  length  (approximately  x/2),  a  is  the  antenna  radius 
and  b  is  the  radius  of  the  plasma  column.  The  other  symbols  and 
functions  are  the  usual. 


THE  INVISCID  FLOW  FIELDS 

Some  estimate  of  the  parameter  top  which  characterizes  the  cold 
plasma  must  now  be  found.  The  simplified  stream-tube  method  of 
Seiff  and  Whiting,2  which  assumes 

(1)  isentropic  flow  along  streamlines,  (2)  flow  parallel  to 
vehicle  axis  and  (3)  the  shape  of  the  pressure  distribution  normal 
to  the  body  is  predictable  (by  first-order  blast  wave  theory).  In 
our  trajectory,  conditions  (1)  and  (3)  are  probably  valid  but  (2)  is 
never  satisfied.  The  angle  of  attack  is  small,  so  we  hesitatingly 
ignore  the  crossflow.  We  also  assume  an  equilibrium  flow  for  all 
altitudes,  which  is  not  exactly  the  case. 

Using  the  profiles  calculated  from  the  inviscid  flow,  we  find 
a  monotonically  increasing  attenuation  from  230  to  100  kft.  (A 
traveling  wave  form  of  current  is  assumed  on  the  antenna.)  Referring 
to  the  flight  data  of  Fig.  5  we  see  this  is  not  the  case. 


THE  VISCOUS  BOUNDARY  LAYER 

Maximum  attenuation  at  193  kft.  indicates  a  shock -boundary 
layer  interaction.  Due  to  the  question  as  to  the  nature  of  the  flow 
and  the  great  boundary  layer  thickness  at  this  altitude,  a  quanti¬ 
tative  calculation  is  subject  to  question.  Nevertheless,  we  proceed 
with  the  phenomenological  approach  to  the  laminar  thermal  boundary 
as  given  by  Schlichting3  since  some  estimate  of  its  effect  is  pref¬ 
erable  to  none. 

The  quantity  that  controls  the  thermal  boundary  layer  overshoot 
is  the  Eckert  number.  We  assume  a  constant  temperature  (300OK)  cyl¬ 
inder.  The  inviscid  flow  calculations  indicate  that  the  kinetic 
temperature  at  the  surface  remains  approximately  constant  at  3000°K 
throughout  the  trajectory.  Since  the  velocity  is  constant  then  the 
Eckert  nimber  is  primarily  dependent  on  the  specific  heat,  which 
for  air  is  a  rapidly  varying  function  of  density  in  the  2500-5000°K 
range.4  To  effect  our  calculations  we  assume  that  the  specific  heat 
of  the  boundary  layer  is  specified  by  the  Inviscid  flow  body  stream¬ 
line.  (A  rather  severe  approximation  due  to  the  variation  of  cp 
through  the  boundary  layer.)  v 
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Using  the  matching  point  scheme  of  Evans  and  Hubert  the  joined 
boundary  layer  -  inviscid  solution  at  140  kft.,  shown  in  Fig.  6,  is 
representative  of  the  results  obtained  at  other  stations  along  the 
cylinder  and  altitudes.  Also  shown  in  rig.  6  is  the  wD  profile  ob¬ 
tained  from  equilibrium  conditions.  w 


DISCUSSION 

From  other  work  based  on  plane  wave  theory  it  was  found  that 
thin  symmetrical  inhomogeneous  layers  could  be  represented  quite 
well  by  equivalent  uniform  layers.  It  was  observed  that  the 
equivalent  thickness  was  well  approximated  by  the  distance  between 
the  two  points  where  w  %  .7  uipmax. 

With  the  v  and  found  from  the  boundary  layer  results  and  the 
equivalent  plasma  radius  determined  as  just  indicated,  Eqs.  (1)  and 
(2)  were  evaluated  for  eight  different  altitude  points.  The  atten¬ 
uation  due  to  the  boundary  layer  is  shewn  in  Fig.  5.  To  effect  better 
agreement  with  experiment,  one  must  add  the  attenuation  due  to  the 
shock  layer  ranging  from  0.2  dB  at  200  kft.  to  4.1  dB  at  120  kft.  A 
more  exact  result  would  be  given  by  a  modification  of  Eq.  (1)  to  in¬ 
clude  two  coaxial  layers. 

In  conclusion,  it  was  found  that  the  VHF  transmission  by  a 
dipole-type  antenna  imbedded  in  an  electrically  thin  axisymmetric 
plasma  sheath  is  explained  by  the  viscous  boundary  layer.  Although 
the  flow  is  quite  complicated,  a  constant  temperature  cylinder  is 
assumed,  and  the  dipole  is  unbalanced,  it  is  believed  that  the  flight 
data  and  analysis  tend  to  show  that  the  antenna  currents  are  somewhat 
isolated  from  the  plasma  by  the  cool  leading  edge  of  the  boundary  layer 
and  the  presence  of  a  dielectric  ablative  coating  over  a  significant 
portion  of  the  antenna.  This  conclusion  is  consistent  with  three  experi¬ 
mental  observations: 

1.  Only  as  the  greatest  attenuation  occurred  did  the  antenna 
resonate. 

2.  There  was  no  significant  pattern  change  as  evidenced  by  the 
ubiquitous  notch. 

3.  The  driving  point  impedance  (Fig.  3)  does  not  indicate 
tight  coupling. 
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ABSTRACT 

The  effect  of  electron  acoustic  wave  generation  has  been  in¬ 
vestigated  for  a  short  cylindrical  dipole  antenna  in  an  aniso¬ 
tropic  (magneto-active)  warm  plasma,  considering  the  collisions 
between  electrons  and  heavy  particles.  Calculations  have  been 
carried  out  for  a  homogeneous  plasma  using  the  hydrodynamic  ap¬ 
proximation  and  the  quasi  static  approximation.  It  is  obtained 
that  the  resistance  due  to  electron  acoustic  wave  increases 
steeply  near  the  upper  hybrid  resonance  frequency.  The  theory 
developed  in  this  paper  can  explain  the  experimental  results 
around  the  above  frequency  range  in  terms  of  emission  of  electron 
acoustic  waves  and  electron  collisions. 


INTRODUCTORY  REMARKS 

The  impedance  of  an  antenna  immersed  in  a  warm  plasma  has 
been  investigated  extensively  in  the  recent  decade.  Considerable 
progress  has  been  made  in  the  case  of  the  warm  isotropic  plasma 
(a  bibliography  will  be  cited  in  a  full  veroxon  of  this  work) .  In 
the  case  of  the  magneto-active  plasma,  however,  the  theory  to 
explain  the  electron  acoustic  wave  effects  on  RF  plasma  probe 
experiments  (Oya  and  Obayashi,  1966,  1967;  Stone  et  al  1966), 
which  has  been  carried  out  in  the  ionospheric  plasma,  is  not 
adequately  developed. 

THEORETICAL  PROCEDURE 


In  this  paper  the  problem  is  investigated  from  the  standpoint 
of  an  RF  probe.  The  short  cylindrical  dipole  antenna  impedance 
is  calculated  in  a  magnetoactive  warm  plasma;  the  EM  field  is 
neglected  since  the  problem  is  restricted  in  a  region  where  the 
electrostatic  field  is  completely  dominant.  The  existence  of' the 
ion  sheath  is  also  neglected  since  the  frequency  range  of  interest 
is  near  the  electron  plasma  and  upper  hybrid  resonance  frequency. 
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The  macroscopic  hydro-dynamic  equation  is  used  to  descrive 
electron  motions  considering  the  electron  collisions  between 
heavy  particles.  Maxwell's  equation,  the  equation  of  continuity 
and  the  equation  of  the  state  are  used  being  transformed  into 
Fourier  space  with  respect  to  the  three  dimensional  problems. 

All  time  functions  are  expressed  by  the  factor  eJ^t  where  m  and 
t  are  the  angular  frequency  and  time  respectively.  In  this  case, 
whole  set  of  basic  equations  is  expressed  in  a  comprehensive  form 
by  introducing  a  dielectric  matrix  which  has  elements  expressed 
as  functions  of  cu,  wave  vector  5,  plasma  frequency  rp.  gyro- 
frequency  fH  and  temperature  T. 

The  basic  equation  is  solved  to  obtain  the  potential  function 
t  as  the  first  step  of  impedance  calculations;  in  this  case  quasi 
static  approximation  is  employed.  One  of  important  boundary 
condition  to  solve  the  differential  equation  is  charge  qs,  and 
current  Js  distributions  on  the  antenna.  It  is  assumed  that 
the  particles  in  the  plasma  are  completely  reflected  back  through 
a  contact  region  between  the  plasma  and  the  RF  probe,  which,  in 
Fourier  transformed  space,  effectively  gives  -jk-Js  +  jwqL  =  0 
on  the  antenna.  Hence, once  we  find  the  charge  function  qg,  the 
problem  can  be  solved  by  obtaining  the  inverse  Fourier  transfor¬ 
mation.  In  this  work,  q~  is  assumed  to  be  distributed  uniformly 
on  a  short  cylindrical  dipole  whose  axis  is  parallel  to  the  DC 
magnetic  field.  This  is  equivalent  to  the  model  of  Balmain 
[1965],  which  is  used  to  obtain  the  impedance  in  the  case  of  an 
isotropic  plasma;  and  the  same  technique  is  .ilso  used  to  obtain 
the  antenna  impedance  by  using  the  e.m.f.  method. 

The  impedance  value  is  finally  expressed  by  an  integration 
in  which  an  integrant  is  a  complicated  complex  function.  This 
integration  is  carried  out  numerically  for  the  ionospheric  plasma 
parameters . 

RESULTS  AND  DISCUSSIONS 


A  numerical  calculation  of  the  impedance  value  is  obtained 
for  a  short  cylindrical  dipole  antenna  whose  elements  have  length 
L  *  1m,  and  radius  pj_  =  0.01m.  In  this  result,  the  gyroresonance , 
the  plasma  resonance,  and  the  upper  hybrid  resonance  are  indicated. 
The  remarkable  evidence  of  the  electron  thermal  effect  is  the 
presence  of  the  acoustic  wave  resistance  which  is  enhanced  near 
the  upper  hybrid  resonance  frequency  and,  also,  in  a  frequency 
range  below  the  gyroresonar  ■*  frequency.  In  a  frequency  range 

higher  than  »,/2  f2  .  *2  the  electron  acoustic  wave  cannot  exist 

pH* 

since  the  wave  length  becomes  smaller  than  the  Debye  shielding 
length;  thus  the  acoustic  resistance  is  not  present. 


In  practice,  we  cannot  neglect  the  existence  of  an  ion 
sheath  surrounding  the  RF  probe,  except  for  the  case  of  sheath 
collaps  condition  which  can  be  artificially  produced  by  applying 
a  high  DC  voltage  to  the  probe.  In  the  frequency  range  in  which 
| Z (a>)  | »  IVCgtt)  i  where  Cs  is  the  so  called  sheath  impedance, 
however,  the  effect  of  sheath  impedance  is  completely  negligible 
[Oya,  1967].  Thus  we  can  confirm  that  above  the  sheath  resonance 
frequency  the  ion  sheath  effect  is  not  serious,  and  around  the 
plasma,  and  upper  hybrid  resonance  frequency  the  ion  sheath  effect 
is  completely  negligible. 

A  microscopic  approach  to  electron  acoustic  waves  in  an 
anisotropic  warm  plasma  was  considered  by  Bernstein  [1958]. 
Bernstein's  relation  is  developed  for  a  case  near  the  upper  hybrid 
resonance  by  Fejer  and  Calvert  [1964].  The  dispersion  equation 
is  only  slightly  different  from  Fejer  and  Calvert's  equation  near 
the  upper  hybrid  resonance  frequency  and  indicate  that  the  theo¬ 
retical  approach  of  this  work  in  a  frequency  range  around  the 
upper  hybrid  resonance  obtains  results  similar  to  the  microscopic 
treatment  of  the  plasma  and  is  sufficient  to  discuss  the  quantita¬ 
tive  nature  of  the  upper  hybrid  resonance. 

In  Fig.  1,  examples  of  impedance  value  using  absolute  value 
and  phase  angle  are  illustrated  for  plasma  parameters  of  the 
ionosphere  at  altitude  of  100  km (a) ,  500  km(b)  and  1000  km(c) . 

The  phase  variation  indicates  that  near  the  upper  hybrid  resonance 
there  is  strong  enhancement  of  the  electron  acoustic  wave;  the 
power  dissipation  of  this  acoustic  waves  affects  the  sharpness 
of  the  resonance.  In  the  collision  dominant  regime  (see  a), 
however,  the  power  dissipation  is  caused  by  the  collisions.  The 
sharpness  of  the  resonance  is  investigated  for  the  ionospheric 
plasma  in  altitude  range  of  100  km  ~  1000  km  being  compared  with 
rocket  observation  result  [Oya,  1967].  There  is  good  agreement 
between  the  theoretical  calculation  and  the  experimental  result. 

CONCLUDING  REMARKS 


The  theory  is  developed  in  an  anisotropic  warm  plasma  for  a 
short  cylindrical  dipole  RF  probe  immersed  in  an  homogeneous 
ionospheric  plasma.  Near  the  plasma  frequency  and  upper  hybrid 
resonance  frequency  the  results  are  completely  useful  and  a  com¬ 
plicated  feature  of  the  experimental  result  obtained  in  the 
ionosphere  can  be  explained.  This  work  may  be  a  starting  point 
for  investigations  of  the  electron  acoustic  wave  effect  in  wider 
frequency  ranges  including  the  cyclotron  harmonic  effects  and 
ion  sheath  effect. 


FREQUENCY  (MHz) 


B  ■  1^  Three  examples  of  calcu¬ 
lated  impedance  value  around 
the  upper  hybrid  resonance 
frequency,  for  the  following 
cases : 


! 

fp(MHz) 

|  fH(MHz) 

T(°K) 

V 

a) 

2.84 

b) 

4.00 

1 

1300 

c) 

1.28 

1  0.96 

1400 
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DETERMINATION  OF  ELF/VLF  TRANSMITTING 
ANTENNA  PERFORMANCE  IN  THE  IONOSPHERE 
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An  experiment  is  planned  that  will  determine  the 
effects  of  the  magnetoionic  medium  on  propagation 
between  an  ELF/VLF  source  in  the  ionosphere  and  a 
surface-based  receiver.  Proper  implementation  of 
this  propagation  experiment  requires  that  instru¬ 
mentation  to  measure  the  driving  point  impedance 
and  efficiency  of  transmitting  electric  dipoles  and 
loop  antennas  in  the  anisotropic  ionosphere. 

An  ELF/VLF  propagation  experiment  is  being  planned  by  the 
Navy  to  determine  how  energy  emitted  from  a  spacecraft  in  the 
ionosphere  propagates  through  the  ionosphere,  couples  across  the 
air-ionosphere  boundary,  energizes  the  Earth-air-ionospheric 
cavity,  and  propagates  to  the  surface  of  the  Earth.  A  high  power 
transmitter  in  an  orDiting  spacecraft  will  be  the  source  of  power. 
Receiving  and  analysis  equipment  at  many  geographic  locations 
will  determine  signal  strength,  polarization  and  direction  of  arrival. 
Total  time  delay  and  frequency  shift  will  be  determined  by  com¬ 
parison  with  telemetered  replicas  of  the  transmitted  ELF/VLF  signals 

However,  losses  and  other  characteristics  of  the  propagation 
path  cannot  be  determined  from  the  received  signals  unless  the 
characteristics  of  the  emitted  signals  are  known.  Some  information 
relative  to  the  radiated  signals  can  be  deduced  from  measurements 
made  on  the  transmitting  antenna.  These  will  be  outlined  in  this 
paper.  In  add:,tion,  the  antenna  measurements  will  provide  a  base 
for  future  experimentation  and  guide  the  design  of  efficient  couplers 
between  the  high  power  transmitting  source  and  the  terminals  of  an 
antenna  immersed  in  the  ionosphere. 

In  order  to  obtain  a  wide  range  of  propagation  path  parameters 
and  many  "in-situ"  conditions  for  measuring  antenna  performance, 
a  polar  orbit  with  an  apogee  of  2500  kilometers  and  a  perigee  of  350 
kilometers  is  planned.  Including  diurnal  variations,  this  will  pro¬ 
vide  a  range  of  electron  densities  from  10^  electrons/meter*  to 
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10  electrons/meter  .  The  predominate  ion  will  change  from  H+  to 
0+  between  the  apogee  and  perigee.  The  other  magnetoionic  para¬ 
meter  affecting  the  antenna  will  be  the  strength  of  the  Earth's  mag¬ 
netic  field  which  will  produce  electron  gyro  frequencies  of  350  kHz 
to  1500  kHz. 


FIGURE  1 .  Angles  determining 
sheath  characteristics 


It  is  anticipated  that  the  ion  sheath  around  the  dipole  elements 
will  be  the  controlling  parameter  affecting  the  antenna  reactance  at 
low  drive  levels.  This  ion  sheath  will  be  perturbed  by  the  vxB 
potentials  induced  in  the  antenna  structure.  The  size  and  shape  of 
the  ion  sheath  will  depend  greatly  upon  the  magnitude  of  v  and  B 
and  the  angles  6  ,  c  ,  and  £ ,  shown  in  Figure  1 ,  which  uniquely 
define  the  orientation  of  the  antenna  relative  to  B  and  v. 

A  recent  analysis  by  Mr.  Dennis  Baker  of  the  Naval  Research 
Laboratory  based  on  such  a  perturbated  sheath  confirms  the  antenna 
capacity  measurements  of  the  LOFTI  II  A  spacecraft  experiment 
which  carried  a  simple  dipole.  The  spinning  antenna  on  LOFTI  II  A 
is  depicted  in  Figure  2  as  wheel  spokes  in  the  x-z  plane.  The  spin 
axis  S  is  along  the  y  axis  and  the  x-y  plane  contains  the  Earth  mag¬ 
netic  vector  3.  The  spacecraft  carried  a  two-axis  magnetometer 
measuring  components  of  the  magnetic  field  parallel  to  the  spin 
axis  and  along  the  antenna.  When  the  antenna  was  in  the  direction 
of  the  x-axis,  the  magnitude  of  the  component  of  the  magnetic  Held 
parallel  to  the  antenna  is  at  a  maximum .  During  the  course  of  the 
experiment  the  angular  displacement  between  the  position  for 
minimum  capacity  and  the  x-vxis  varied.  Among  other  results  of 
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ABSTRACT 

The  near-field  phase  and  amplitude  broadside  to  a  long  hori¬ 
zontal  dipole  near  Byrd  Station,  Antarctica,  are  measured  at  12.8 
khz.  The  resulting  data,  together  with  an  appropriate  mathemati¬ 
cal  model  suggest  reasonable  values  of  depth  and  bulk  complex 
permittivity  for  the  glacial  ice,  and  a  conductivity  figure  for 
the  underlying  earth. 

INTRODUCTION 

The  phase  and  amplitude  characteristics  of  the  near  field 
broadside  to  a  long  horizontal  dipole  at  very  low  frequencies 
constitute  a  potential  source  of  useful  information  about  the 
nature  of  the  underlying  media  [Biggs,  1968].  This  study  concerns 
Antarctic  glacial  terrain,  and  makes  use  of  data  collected  re¬ 
cently  in  the  Byrd  Glacier  area  to  draw  tentative  conclusions 
about  the  nature  of  the  sub-surface  media.  The  method  suggests 
the  depth  of  the  ice  shelf  is  in  the  neighborhood  of  2250  meters, 
that  the  bulk  relative  complex  dielectric  constant  of  the  ice  may 
be  approximated  by  the  figure  6.0  -  jl4,  and  that  the  conductivity 
of  the  underlying  earth  approaches  0.10  mho  per  meters.  All  of 
these  figures  are  shown  to  be  reasonable  in  light  of  available 
data  from  other  sources. 

EXPERIMENT  AND  MATHEMATICAL  MODEL 

Two  electric  dipoles  21  and  8  miles  long  respectively,  lie  at 
right  angles  near  the  top  of  a  thick  ice  shelf  at  an  Antarctic 
site  12  miles  from  Byrd  Station;  they  may  be  operated  indepen¬ 
dently  of  each  other  in  the  3-30  khz.  frequency  range.  Although 
they  are  used  primarily  for  polar  ionospheric  studies,  the  follow¬ 
ing  auxiliary  experiment  was  conducted  during  January  1969.  With 
the  21  mile  dipole  excited  at  12.8  khz.  and  the  8  mile  dipole  un¬ 
excited,  a  vehicle  equipped  with  VLF  receiving  and  telemetry 
instrumentation  was  driven  away  from  the  dipole  center  along  a 
broadside  path  with  respect  to  the  axis  of  the  21  mile  antenna. 

It  sampled  the  vertical  magnetic  field  over  a  total  distance  of 
11  miles,  approximately  3/4  wavelength,  in  the  general  direction 
of  Byrd  Station.  The  field  amplitude  was  measured  along  this  path 
at  ground  level;  simultaneously,  the  signal  was  telemetered  back 
to  the  station  for  phase  comparison  at  the  transmitter.  These 
data  are  interpreted  in  terms  of  theoretical  predictions  to  be 
described.  Further  comparison  is  made  with  findings  I’ecently 
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will  be  under  control  to  provide  for  variation  from  low  levels  up  to 
the  power  obtainable  from  a  source  capable  of  supplying  better  than 
1  kilowatt  into  the  antenna  coupling  unit.  This  power  will  be  mo¬ 
dulated  to  provide  a  constant  envelope  signal  or  various  pulse  com¬ 
binations  .  This  provides  the  experimenter  with  a  wide  range  of 
parameters  for  propagation  path  experimentation,  ionospheric 
excitation,  and  antenna  performance  determination. 

The  spacecraft  will  contain  a  broad  spectrum  of  instrumentation 
to  measure  .he  parameters  affecting  the  experiment;  in  particular, 
methods  of  measuring  absolute  spacecraft  attitude,  antenna  orient¬ 
ation  relative  to  the  Earth's  magnetic  field,  and  devices  for 
measuring  electron  density  and  ion  constituency. 

A  prime  factor  to  be  known  about  antennas  is  the  far-field 
pattern.  For  this  purpose,  the  ground  receiving  system  is  quite 
limited.  Actually,  any  practical  measurement  technique  will  be  a 
compromise.  However,  other  experimenters  will  have  orbiting 
spacecraft  containing  receiving  equipment  which  will  provide  some 
indication  as  to  the  far  field  In  the  ioncsphore.  Similarly,  one  or 
more  spacecraft  launched  and  separated  from  the  transmitting  space¬ 
craft  could  provide  information.  A  receiving  rocket  launched  in  a 
trajectory  cutting  through  the  radiated  field  would  also  provide  field 
pattern  information.  In  addition,  a  receiver  in  the  spacecrat.'-  acting 
lUe  a  top-side  sounder  would  establish  confidence  that  the  signal 
vas  radiated  by  receiving  energy  reflected  or  refracted  back  from 
discontinuities  or  inhomogenuities .  Unfortunately,  it  is  impractical 
to  thoroughly  scan  the  entire  volume  to  determine  the  patterns ,  how 
they  change  with  radial  distances,  and  the  relative  power  densities 
in  the  several  propagating  modes. 

Most  information  relative  to  antenna  performance  will  be  obtained 
at  the  transmitting  spacecraft. 


FIGURE  3 . 
Impedance 
measurement 
points 
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First,  the  current  into  the  antenna  measured  at  B  (shown  in 
Figure  3)  and  the  voltage  measured  at  A  across  the  antenna  termi¬ 
nals,  together  with  their  relative  phases,  will  be  telemetered  to 
Earth.  This  data  will  be  correlated  with  the  environmental  data  and 
experiment  parameters. 

Secondly  (as  shown  in  Figure  3),  the  antenna  coupling  unit  and 
matching  transformer  will  be  servo-controlled  to  obtain  maximum 
power  transfer.  The  position  of  the  transformer  tap  and  the  reactance 
position  will  be  telemetered  to  monitor  the  performance  of  the  servo 
system.  This  information,  providing  the  range  and  bandwidth  cf  the 
servo  system  is  adequate  to  follow  the  changes ,  will  confirm  the 
results  obtained  by  measurement  of  voltage,  current,  and  phase. 

It  is  anticipated  that  there  will  be  non-linear  effects  observ¬ 
able  by  changing  the  level  of  the  drive  power.  Harmonics  and  other 
frequencies,  detectable  on  the  ground  or  in  the  spacecraft  receiver, 
will  be  generated  and  related  to  the  antenna  performance  behavior. 
This  non-linearity  will  be  reflected  in  the  waveforms  of  the  antenna 
current  and  driving  voltage.  These  waveforms  will  likely  be  tele¬ 
metered  to  Earth  but  only  for  special  tests  since  they  require  large 
bandwidths  in  the  telemetry  channels.  In  order  to  conserve  tele¬ 
metry  bandwidth,  these  waveforms  (as  shown  in  Figure  4)  will  be 
processed  by  a  harmonic  analyzer  on  the  spacecraft  and  the  results 
relayed  in  digital  form  to  the  telemetry  stations . 
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FIGURE  4 . 
Method  of 
determining 
non-linear 
characteristics 


Potentials  will  probably  be  sufficiently  large  to  cause  the  sheath 
to  collapse  at  one  end  of  the  antenna  as  shown  in  Figure  5.  Simple 
measurement  of  reactance  across  the  antenna  terminals  A  and  B  will 
not  be  sufficient  to  reveal  this  asymmetry.  Measurement  of  the 
capacity  between  each  of  the  antenna  terminals  and  the  body  of  the 
spacecraft,  terminal  C,  will  provide  some  information  relating  the 
capacity  between  the  sheath  and  each  of  the  dipcle  elements . 
Additional  information  relative  to  the  sheath  formation  will  be  ob¬ 
tained  by  putting  direct  current  biases  between  the  antenna  elements 
and  the  spacecraft  and  by  an  electron-gun  on  the  spacecraft. 
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FIGURE  5 . 
Measurement  of 
pulsating  sheath 


Figure  5  shows  a  static  representation  of  the  sheath.  When  the 
antenna  is  used  for  transmission,  the  driving  potential  across  the 
antenna  terminals  and  along  the  antenna  may  equal  or  exceed  the 
v  x  B  potentials.  In  any  case,  there  will  be  a  sheath  which  oscil¬ 
lates  from  one  end  of  the  dipole  to  the  opposite  end  of  the  dipole  at 
the  driving  frequency.  Measurements  at  the  driving  frequency  will 
not  reveal  the  full  dynamic  nature  of  the  sheath.  In  order  to  study 
this  pulsating  sheath,  these  capacities  will  be  measured  at  fre¬ 
quencies  of  the  order  of  ten  to  a  hundred  times  that  of  the  driving 
frequency.  The  representation  of  the  circuit  by  three  capacitors  in 
the  lower  part  of  Figure  5  is  grossly  simplified,  since  the  capacity 
is  distributed  and  the  ionosphere  bounding  the  sheath  is  not  an 
electrical  short  but  has  a  complex  impedance  as  diagrammatically 
indicated  in  the  upper  part  of  the  figure.  In  order  to  determine  some 
characteristics  of  the  antenna-sheath-plasma  circuit,  these  im¬ 
pedances  will  be  measured  at  several  frequencies  and  a  sheath 
impedance  derived  from  the  data . 

This  experiment  is  still  in  the  formulation  stages  and  should 
be  of  great  interest  to  those  interested  in  the  effect  of  plasma  and 
anisotropic  media  on  antennas.  Any  comments  to  make  this  a  more 
definitive  and  productive  experiment  would  be  appreciated. 
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ABSTRACT 

The  near-field  phase  and  amplitude  broadside  to  a  long  hori¬ 
zontal  dipole  near  Byrd  Station,  Antarctica,  are  measured  at  12.8 
khz.  The  resulting  data,  together  with  an  appropriate  mathemati¬ 
cal  model  suggest  reasonable  values  of  depth  and  bulk  complex 
permittivity  for  the  glacial  ice,  and  a  conductivity  figure  for 
the  underlying  earth. 

INTRODUCTION 

The  phase  and  amplitude  characteristics  of  the  near  field 
broadside  to  a  long  horizontal  dipole  at  very  low  frequencies 
constitute  a  potential  source  of  useful  information  about  the 
nature  of  the  underlying  media  [Biggs,  1968].  This  study  concerns 
Antarctic  glacial  terrain,  and  makes  use  of  data  collected  re¬ 
cently  in  the  Byrd  Glacier  area  to  draw  tentative  conclusions 
about  the  nature  of  the  sub-surface  media.  The  method  suggests 
the  depth  of  the  ice  shelf  is  in  the  neighborhood  of  2250  meters, 
that  the  bulk  relative  complex  dielectric  constant  of  the  ice  may 
be  approximated  by  the  figure  6.0  -  j!4,  and  that  the  conductivity 
of  the  underlying  earth  approaches  0.10  mho  per  meters.  All  of 
these  figures  are  shown  to  be  reasonable  in  light  of  available 
data  from  other  sources. 

EXPERIMENT  AND  MATHEMATICAL  MODEL 

Two  electric  dipoles  21  and  3  miles  long  respectively,  lie  at 
right  angles  near  the  top  of  a  thick  ice  shelf  at  an  Antarctic 
site  12  miles  from  Byrd  Station;  they  may  be  operated  indepen¬ 
dently  of  each  other  in  the  3-30  khz.  frequency  range.  Although 
they  are  used  primarily  for  polar  ionospheric  studies,  the  follow¬ 
ing  auxiliary  experiment  was  conducted  during  January  1969.  With 
the  21  mile  dipole  excited  at  12.8  khz.  and  the  8  mile  dipole  un¬ 
excited,  a  vehicle  equipped  with  VLF  receiving  and  telemetry 
instrumentation  was  driven  away  from  the  dipole  center  along  a 
broadside  path  with  respect  to  the  axis  of  the  21  mile  antenna. 

It  sampled  the  vertical  magnetic  field  over  a  total  distance  of 
11  miles,  approximately  3/4  wavelength,  in  the  general  direction 
of  Byrd  Station.  The  field  amplitude  was  measured  along  this  path 
at  ground  level;  simultaneously,  the  signal  was  telemetered  back 
to  the  station  for  phase  comparison  at  the  transmitter.  These 
data  are  interpreted  in  terms  of  theoretical  predictions  to  be 
described.  Further  comparison  is  made  with  findings  recently 
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obtained  by  the  U.S.  Army  Cold  Regions  Research  Laboratories  re¬ 
garding  the  underlying  ice  structure  in  connection  with  a  six  inch 
core  drilled  through  the  glacier  at  Byrd. 

The  relevant  mathematical  model  is  illustrated  in  Figure  1; 
it  comprises  a  homogeneous,  isotropic,  lossy  dielectric  ice  layer 
covering  a  semi-infinite  region  of  imperfectly  conducting  earth. 
The  ice  layer  lies  below  a  semi-infinite  region  whose  properties 
are  those  of  free  space;  plane  boundaries  are  assumed  throughout. 
Assuming  a  horizontal  line  source  at  the  air-ice  interface,  ex¬ 
pressions  may  be  derived  for  the  horizontal  electric  field  at  this 
boundary  in  terms  of  distance  from  the  source.  These  are  well 
known,  and  have  been  presented  before  [Wait,  1962;  Collin,  I960; 
Biggs  and  Swarm,  1965].  The  vertical  magnetic  field  is  of  course 
simply  related  to  these  expressions  through  Maxwell’s  equations. 
Digital  computer  techniques  are  used  to  obtain  accurate  solutions; 
they  are  cast  in  a  general  form  to  permit  application  to  a  variety 
of  natural  media,  of  which  ice  and  underlying  earth  represent  a 
special  case. 


21  mile  dipole 


oir  region  (free  space) 


ice  layer 
^lce»H°) 


Figure  f. 
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It  is  shown  that  two  basic  modes  are  the  primary  contributors 
to  “he  total  field  at  the  interface,  namely  a  component  represent¬ 
ing  propagation  in  the  ice  layer  and  one  propagating  along  the 
interface  itself.  An  additional  weak  lateral  wave  contribution 
results  from  propagation  just  below  the  earth-ice  interface.  A 
brief  mathematical  discussion  summarizes  their  origins  in  terms 
of  an  integral  evaluation  in  the  complex  plane.  Predictions  are 
made  regarding  the  phase  and  amplitude  of  the  electric  field  at 
the  air-ice  interface  in  terms  of  a  range  of  possible  complex 
dielectric  constants  and  layer  depths.  The  best  fit  to  the 
experimental  data  is  found  to  occur  for  the  combination  of  param¬ 
eters  indicated  above.  Figures  2  and  3  illustrate  the  agreement 
obtained  under  these  conditions. 
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DISCUSSION  AND  CONCLUSIONS 

The  vertical  thickness  of  the  glacier  at  Byrd  Station  has 
been  measured  directly  by  drilling  to  the  bottom  [Gow,  Ueda  and 
Garfield,  1968];  the  resulting  depth  was  2164  meters,  which  com¬ 
pares  well  with  our  value  of  2250  meters.  Other  investigations 
have  yielded  depth  figures  varying  up  to  2570  meters  in  the 
neighborhood  of  Old  Byrd  Station  [Thiel,  Bentley,  Ostenso,  and 
Behrendt;  1959],  where  the  ice-rock  interface  is  believed  to  be 
quite  irregular  along  a  path  roughly  parallel  to,  and  only  a  few 
miles  from  the  path  traveled  by  our  vehicle. 

A  conductivity  of  0.1  mho  per  meter  is  slightly  high  for  the 
underlying  rock,  but  not  unreasonable  for  the  mixture  of  granite, 
ice  and  water  found  at  the  bottom  of  the  drill  hole  by  Gow,  Ueda 
and  Garfield.  Furthermore,  our  method  is  relatively  insensitive 
to  the  assumed  conductivity  of  the  underlying  medium,  so  long  as 
this  is  greater  than  0.001  mho  per  meter. 

No  experimental  data  are  available  regarding  the  complex 
permittivity  of  the  ice  as  a  function  of  depth  at  the  drill  site 
or  elsewhere  in  the  Byrd  Glacier  area,  although  the  density  and 
complex  dielectric  constant  of  the  surface  snow  have  been  measured 
at  several  different  temperatures  at  10  khz.  [Rogers,  1967];  e  = 
eQ(2.7  -  j2.7)  at  -  30°  C,  with  a  density  of  approximately  0.4 

gm/cc  is  a  typical  figure;  the  complex  permittivity  decreases  with 
temperature ,  and  increases  with  pressure  packing  until  impermeable 
ice  results  at  a  depth  of  approximately  65  meters  [Hatherton, 
1965],  Assuming  that  this  parameter  is  essentially  independent  of 
pressure  over  the  entire  vertical  depth  of  the  glacier,  and  pri¬ 
marily  dependent  on  temperature  and  the  presence  of  impurities,  it 
is  possible  to  infer  some  bounds  on  the  expected  average  figure. 
Simple  straight-line  approximations  to  the  temperature  profile 
measured  in  the  drill  hole  [Gow,  et.  al.]  yield  an  average  tem¬ 
perature  near  -20°  C.  This  temperature,  in  turn,  suggests  a  lower 
bound  of  3  -  j6,7  for  the  complex  permittivity  in  terms  of  the 
figure  for  pure  ice  at  12.8  khz.  at  this  temperature  [Watt  and 
Maxwell,  I960],  Impurities  in  the  glacial  ice  can  be  expected  to 
increase  the  dielectric  and  loss  factors.  A  s»*ep  toward  account¬ 
ing  for  their  presence  is  made  by  assuming  them  to  be  comparable 
to  or  less  than  those  of  the  Athabasca  Glacier  in  Alberta,  Canada, 
for  which  Watt  and  Maxwell  give  some  data.  Their  figure  of  14 
-j21,  measured  in  the  surface  ice  at  0°  C,  is  taken  as  an  upper 
bound  since  it  is  known  that  the  quantity  of  interest  decreases 
at  lower  temperatures.  Moreover,  the  greater  isolation  of  the 
Byrd  Glacier  from  civilization  suggests  a  lower  impurity  level 
relative  to  that  at  the  Athabasca  Glacier.  On  these  grounds,  our 
value  of  6  -  jl4  (er  =  6,  tan  6  =  2.33)  appears  justified. 


It  is  concluded  that  the  measurement  of  the  phase  and  ampli¬ 
tude  of  the  VLF  ground  wave  can  provide  information  about  the 
nature  of  glacial  ice  and  the  sub-surface  land  boundary.  It  may 
be  viewed  as  a  valuable  adjunct  to  other  methods  currently  in  use 
to  study  these  characteristics  of  the  Antarctic  terrain. 

Acknowledgement :  The  authors  wish  to  thank  Mr.  J.  C.  Rogers  and 
Mr.  A.  Chandler  of  the  University  of  Washington  Electrical  Engi¬ 
neering  Department  for  helpful  discussions  and  for  assistance  with 
the  in-site  measurements  at  the  U.W.  Longwire  Antenna  site.  This 
work  was  supported  by  the  National  Science  Foundation  under  grant 
No.  GA-161Q. 


REFERENCES 

Biggs,  A. W.  (May  1968),  Geophysical  Exploration  in  polar  areas 
with  very  low  frequency  phase  variations,  IEEE  Trans.  Ant.  6  Prop. 
Vol.  AP-16,  No.  3,  364-365.  "  ' 

Wait,  J.R.  (1962),  Electromagnetic  Waves  in  Stratified  Media, 
MacMillan  Book  Company,  New  York. 

Collin,  R.E.  (1960),  Field  Theory  of  Guided  Waves,  McGraw-Hill 
Book  Company,  New  York. 

Biggs,  A.W.  and  Swarm,  H.M.  (1965),  Analytical  Study  of  the 
Radiation  Fields  from  an  electric  dipole  in  stratified  and  in¬ 
homogeneous  Antarctic  terrain.  Technical  Report  No.  98,  Dept,  of 
E.E.,  Univ.  of  Wn. 

Gow,  A.J.,  Ueda,  H.T. ,  and  Garfield,  D.E.  (6  Sept.  1968),  Antarctic 
ice  sheet:  preliminary  results  of  first  core  hole  to  bedrock, 

Sci.  161,  1011-1013. 

Rogers,  J.C.,  private  communication 

Hatherton,  T.  (1965),  .intarctica,  A.H.  and  A.W.  Reed  Publishers  of 
New  Zealand  Books,  Wellington,  Auckland  N.Z. 

Thiel,  E. ,  Bentley,  C.R.,  Ostenso,  N.A.,  and  Behrendt,  J.C., 

(July  1959),  Oversnow  traverse  programs,  Byrd  and  Ellsworth 
Stations,  Antarctica,  1957-58;  seismology,  gravity,  and  magnetism, 
IGY  Glac.  Rpt.  Series  No.  2,  1GY  World  Data  Center  A,  Glaciology, 
Am.  Geog.  Soc.  Sec.  II,  5-e. 

Watt,  A.D.  and  Maxwell,  E.L.,  (1960),  Measured  electrical  proper¬ 
ties  of  snow  and  glacial  ice,  J.  Res.  NBS  64D  (Radio  Prop.)  No.  4, 
357-63. 

Tai,  C.  T.  (1951),  The  effect  of  a  grounded  slab  on  radiation  from 
a  line  source,  J.  Appl,  Phys.,  22,  405. 


-55- 


RADIATION  FROM  AN  ELEVATED  VERTICAL  DIPOLE 
ABOVE  AN  EXTENDED  GROUND  SCREEN 

W.  T.  PATTON 


RCA,  MOORESTOWN,  NEW  JERSEY 

A  good  inductive  ground  screen  used  to  provide  low  angle 
space  wave  radiation  from  a  vertical  dipole  can  support 
the  Sommerfeld  surface  wave  since  the  phase  angle  of  the 
effective  surface  impedance  can  exceed  45  degrees  and  may 
approach  90  degrees  for  a  very  good  screen.  The  Sommerfeld 
attenuation  functior  is  important  to  accurate  calculation 
of  ground  screen  effects  for  screens  of  large  extent. 


The  ability  of  an  antenna  system  to  establish  long  range  links  via  the 
ionosphere  is  determined  by  the  characteristics  of  the  space  wave 
radiation  pattern  at  low  elevation  angles.  For  a  short  vertical  dipole 
over  a  homogeneous  earth,  the  pattern  is  given  by 


Jkhsin^  .  „  -jkh  sin'J<  .  .  Z,,,. 

e  +  Rv  e  J  sir#  -  1QJQ 

G1  =  cosif)  - - - where  Rv  =  fr 

Z  ,(*) 

simp  — 

is  the  Fresnel  reflection  coefficient  for  vertical  polarization. 

Zi<*>  i  /  2,\in  /  v1/9 


i  L  cos2<iA 1/2  /  o  Xi 

y-  -  —  *  nier  -tj 


(1) 

(2) 

z1(*) 

~¥~ 

(3) 


is  the  surface  impedance  of  a  finitely  conducting  earth.  The  space 
wave  radiation  pattern  (1)  of  a  vertical  dipole  generally  shows  sharp 
attenuation  for  elevation  angles  below  10  to  15  degrees. 


Two  methods  of  improving  low  angle  radiation  are  1)  increasing  vertical 
directivity  by  elevating  the  dipole  and  2)  providing  a  ground  screen 
under  the  antenna  of  sufficiently  large  extent  to  be  effective  at  low 
elevation  angles.  Ground  screens  longer  than  100  wavelengths  are 
required  to  control  the  pattern  at  angles  about  5  degrees  elevation. 
Following  Wait  (1967)  with  the  Sommerfeld  attenuation  function  included 
as  in  Wait  (1963)  the  space  wave  radiation  pattern  of  the  composite 
system  is  giveu  approximately  by 

(4) 

J>  cost))  )  W  (X,Z2)dX 
The  Sommerfeld  attenuation  function  is  (5) 

W(X,Z2)  =  1  -j.J’rrp”  e  W  erfc  (j ^/w),w  =/l  +f  p  ,p  =  “j 
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and  Z 2  is  the  composite  surface  impedance  formed  by  the  parallel 
combination  of  the  surface  impedance  of  the  earth  (equation  3)  and  the 
impedance  of  the  ground  screen.  If  the  ground  screen  is  composed  of 
parallel  wires  or  a  wire  grid,  its  impedance  is  given  approximately  by 


where  d  is  the  grid  spacing  and  c  is  the  wire  radius.  The  grid 
impedance  is  inductive  and,  to  be  effective  in  improving  low  angle 
radiation,  must  be  much  smaller  than  surface  impedance  of  the  ground. 
Thus  the  composite  surface  impedance  determining  the  behavior  of  the 
Sommerfeld  attenuation  function  will  be  inductive,  with  a  phase  angle 
approaching  90  degrees.  The  magnitude  and  phase  of  this  function  shown 
in  Figures  1  and  2  show  the  effect  of  the  Sommerfeld  surface  wave  as 
noticed  by  Wait  (1957)  for  highly  inductive  ground  systems  (phase  angle 
over  45  degrees).  It  is  manifest  by  the  attenuation  function  exceeding 
unit  magnitude  and  undergoing  rapid  phase  variation. 

Satisfactory  results  can  be  obtained  from  equation  (4)  by  setting  W 
equal  to  one  only  when  the  magnitude  of  the  numeric  distance  P  in 
equation  (5)  is  much  less  than  one.  This  may  often  be  the  case  for 
antenna  systems  where  good  radiation  below  elevation  angles  of  about 
10  to  15  degrees  is^f  primary  concern.  However,  for  good  performance 
at  lower  elevation  angles,  an  extended  ground  screen  must  be  used  and, 
in  the  interest  of  economy,  the  largest  wire  spacing  consistent  with 
good  performance  is  desirable.  In  this  case,  the  Sommerfeld  attenuation^ 
function  must  be  included  in  equation  (4). 


The  effect  of  grid  spacing  on  the  elevation  pattern  can  be  determined 
initially  by  considering  the  composite  ground  system  to  be  of  infinite 
extent.  The  patterns  shown  in  Figure  3  are  for  this  case.  The 
interesting  result  is  that  for  large  grid  spacings,  the  radiation 
pattern  of  the  elevated  dipole  is  poorer  than  that  over  unaided  ground. 
The  reason  for  this  can  be  seen  when  equation  (1)  is  approximated  for 
very  small  elevation  angles 


cosiJj  sini]; 

7 

sinijr*-  2 
the  ii^>e 


( 


1  +  jkh 


(7) 


where  the  subscript  2  denotes 

by  the  wire  grid.  Since  is  highly  inductive,  ^ 
dipole  height  h  is  increased  and  obtains  a  minimum  at 


ipedance  of  the  earth  as  augmented 
G„  is  reduced  as  the 

a  height  depending 
upon  the  magnitude  and  phase  of  Z_.  The  pattern  of  a  finite  extended 
ground  screen  will  tend  to  be  better  than  those  shown  in  Figure  3  at  the 
lower  elevation  angles  due  to  diffraction  from  the  end  of  the  screen. 


The  Sommerfeld  attenuation  constant  should  be  included  in  calculations 
leading  to  the  design  and  evaluation  of  economical  ground  systems  for 
low  angle  radiation.  Elevating  the  dipole  above  the  ground  screen  can 
in  some  cases  result  in  poorer  performance  than  might  be  had  with  a 
monopole  on  the  ground. 
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Magnitude  of  Sommerfeld’s 
Attenuation  Function 

Figure  i 


Radiation  From  An  Elevated  Dipole 
Above  A  Composite  Ground  System 


Phase  of  Soramerfeld's 
Attenuation  Function 


Figure  3 


Figure  2 
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NEAR  FIE  ID  SOLUTIONS  FOR  A  VERTIGAL  HERTZIAN  DIPOLE 
OVER  A  FINITELY  CONDUCTING  EARTH 

David  C.  Chang,  University  of  Colorado 
James  R.  Wait,  ESSA  Research  Laboratories,  Boulder,  Colorado  §0302 


Various  approximations  to  the  exact  formal  integral 
representation  are  considered.  A  new  approximate  form 
is  obtained  which  appears  to  be  valid  even  when  the  dipole 
is  near  a  poorly  conducting  earth.  Numerical  compari¬ 
sons  with  the  exact  integral  formula  verify  some  of  the 
conjectures. 

1.  INTRODUCTION 

The  complete  electromagnetic  field  of  a  vertical  Hertzian 
electric  dipole  of  current  moment  IQ  ds,  located  at  a  distance  z0 
above  a  homogeneous  half-space  of  conductivity  a  and  dielectric 
constant  e  can  be  determined  from  a  Hertz  vector  II  which  has 
only  a  z  component  [Banos,  1966].  The  formal  exact  solution  for 
the  time-harmonic  problem  is  expressible  in  the  form 


E-VxVxUa  ;  H  =  (o-i«e)VxII  a  ; 

z  z  z  z 

nz(p,z)  =  (-  V°-  )k(z,  zo;p);  K  =  K(p)  +  K(#)  , 
where  the  primary  influence  is 


(z,  ZQ;p) 


ik  [(z-z  f  +  fl* 
o  o 

a 

k  [<z-z  )2  fr  p2]* 
o  o 


(1) 

(2) 

(3) 


and  the  secondary  influence  is 

.  r  ik  (z+z  )[l-o2] 

K1  *(z,  z  ;  p)  =  i  \  R(«)  a 


2li 


„  , ,  ,  c  d  a 

°  c  [1-c2]* 


(4) 


R(«)  =  [a2  (1-c2 )*  -  (n2  -  a2  } fyfa2  (1-a2 )*  +  (n2  -  a2 )^  ]  .  (5) 

Here,  n  =  [(o  -  i  a)  e)/(-ia>  e  )J^  is  the  (comolexl  refractive  index 
of  the  earth;  £  =  120  tr,  the  characteristic  impedance  of  the  free- 
space;  k  ,  the*  free -space  wave  number  and  aE  is  the  unit  vector 
along  the0  z  direction.  In  (l)-{5),  the  time -dependence  e"iwt  is 


as  sun  ed.  Since  the  secondary  influence  K'  ,  in  (4),  is  expressed 
onlv  in  terms  of  an  integration  in  the  Fourier -transform  plane  a  , 
various  asymptotic  solutions  have  been  developed  for  different  ranges 
of  observation  (Banos,  1966].  Specifically,  the  knowledge  of  the  field 
near  the  vertical  axis  (i.  e. ,  p  =*•  0)  is  important  in  the  determination 
of  the  self -impedance  of  the  source  dipole  using  the  e.  m.  f.  method 
(Wait,  1962]  as  well  as  in  the  finding  of  a  better  approximate  solution 
for  the  current  of  a  finite  linear  antenna  using  the  integral  equation 
approach  (Chang  and  Wait,  1969].  In  the  past,  appr  .  ximate  formulas 
which  give  the  expression  in  closed  forms  have  been  developed  for 
cases  when  the  source  dipole  is  located  far  away  from  the  earth 
(kQ  (z  *  z  j  £  10)  (Jasik,  1961] ,  and  when  the  earth  is  highly  conductive 
[Wait,  1962;  Lavrov  and  Knyazev,  1965],  However,  h  a  more  general 
situation,  the  source  and  the  observation  point  are  near  the  ground 
which,  itself,  may  not  be  highly  conductive.  In  this  case,  the  avail¬ 
able  formulas  fail  to  give  an  adequate  description  of  the  physical 
problem.  Motivation  for  a  more  general  analysis  has  arisen  in 
connection  with  the  design  of  HF  antennas  for  special  applications  and 
more  recently  in  the  remote  sensing  of  the  moon's  surface. 

2.  APPROXIMATE  SOLUTION  NEAR  THE  VERTICAL  AXIS 

In  the  following  derivation,  we  make  the  assumptions: 
p  « (z  +  zQ)  and  kQp«  1  ,  (6);  |n|2  u  10,  (7) 

Branch  cuts  at  a  =  *  1  and  c  =  ±nfl  in  the  expressions  (4)  a^d 
(5)  are  so  chosen  that  jarg  (a2  -  n2)? \s  ir/2  and  jarg  (a2  -  1)*  J  <  n/Z. 
We  now  let  z  +  zQ  =  d  and  split  the  integration  path  in  (4)  into  two 

P***'-  x  . 


,,(s).  .  ./■  P  .  f  «  ,  »  ik_d(l  -  a2)3  T  ,  eda 

K'  (z,z  ;p)  =  H  \  +  \  ;R(a)e  o  J  (k  p  a) - j- 

°  °  °  (l-a2)* 

A  i  J,  A 

=  i  \  [n2a  -  {n2  -  1  +  e2)2]/[n2a  -  {n2  - 1  +  a2)2]  J^pfl  -  «2J  2 

®  CO 

x  eik°^ad«  +  ^*[in2  s  -  (n2  - 1  -  u2)2  ]/[i  n2  c  +  (n2  -  1  -  e2)2  ] 


„ikoda 


da  + 


x  elkodc  JQ(kop[l  4  a2)2)  da 

i 


In  deriving  (8),  a  change  of  variables  (1  -  i2 ) 2  —  a  and  (a2  -  l)2-»  a 
has  been  applied  to  the  first  and  second  integrals,  respectively.  ^ 
In  the  first  integration,  a  s  1,  we  can  approximate  (n2  -  1  +  a2)2 
by  n,  when  it  is  compared  with  n2  c  ,  based  on  (7).  This 
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approximation  apparently  is  better  than  o^jnj  j  .  On  the  other 
hand,  the  second  integrand  contributes  significantly  only  when 
a  £  l/(k  d)  because  of  the  exponential  decay  term  in  the  integrand. 
Again,  based  upon  (7),  we  can  approximate  (n2  -  1  -  c2)fc  by  a  . 
Also,  because  of  (6),  the  Bessel  function  J  (k  pfl  +  a2]^)  can  be 
approximated  by  1.  To  determine  the  error  involved,  we  can  esti¬ 
mate  the  remainder  or  "error"  term: 


“  4 

n.  „  C  r*i  n2  a  -  (n2  - 1  -  azy  i  n  a  - 1 1 

a,"<”  = 

o 

lA>d  A 

a.  C  r  i n2  c  - (n2  - 1  -  a1)  ino-1  ~j 

5  +  ’  in°  +  1J 


k  ad  i 

°  J  (k  p[l  +  o2]*)do 
o  o  1 


da 


It  ie  not  difficult  to  show  that  this  is  of  the  order  of  l/\n  |  "  Ol* 
l/(Jn|*  k2  d2) .  Consequently,  we  can  approximate  (8)  as 

K*8*(z,  zo;  p )  ^  Kc  (z,  zq) 


pi  i  a  k  d 

=s  i  \  [(na-l^(na  +  I)]  e  °  da  +  \  [(i n  c  -l)/(in  a  +  i)]e 


-  a  k  c 
o 


This  can  be  recombined  to  give 

*  ,  ik  d(l-c2)* 

K  fz,  z  )  =  i  \  (1  -  Z/tn{l  -  a2)*  +  I])  e  °  - 1 

c  °  -  (1  -a2)* 


=  K  (z,  z  )  +  K  (z,  z  ) 

Cj  O  Cz  o 


ik  (z+z  1 
o  o 

b 

k  (z  +  z  ) 
o  o 


where  m  i 

P  ik  d(l-c2r  - 

K  (z,  z  )  =  i  \  e  °  — ■ — “1  = 

c:  °  -  (1 -a2)* 

a»d  w  1 

p  1  ik  d(l-c2)* 
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-ik  d/n  ____  _ 

=  — -  e  °  E,  (-ik  r  }  ;  r  i  (z  +  z  )  (1  +  l/n)  .  (11) 

n  o  o 

where  E*  is  the  exponential  integral  of  the  first  order  [Abramowitz 
and  Stegun,  1964] ,  Substitution  of  (10)  and  (11)  into  (9)  yields 

-ik  d/n 

♦  ~e  °  Ei(-ik  r)  ,  (12) 

n  o 


K  (z.z  )  = 
c  o 


ik  (z+z  ) 
o  o 
e 

k  (z  +  z  ) 
o  o 


where  the  complex  distance  r  is  defined  in  (11).  It  is  interesting  to 
note  that  when  jnj  »  1  and  k^d/Jnj  «  l,  (12)  becomes 


K 
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(z,  z  ) 
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2m 


ik  (z+z  ) 
o  o 
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k  (z  +  z  ) 

G  O 


+  —  iM-tk  !*+*,]> 


K.(z.z)  (13) 

■L.  O 


which  agrees  with  the  approximate  solution  obtained  by  Wait  [1962] 
and  Lavrov  and  Knyazev  [19651. 


In  the  limiting  case  when  jnj  —  »  ,  a  perfectly  conducting 
ground,  only  the  first  term  in  (12)  remains  which  corresponds  ex¬ 
actly  to  the  contribution  from  an  image  source  below  the  ground. 
On  the  other  hand,  if  k  (z  +  z  )  »  1  ,  using  the  asymptotic  ex¬ 
pansion  [Abramowitz  ancl  Stegun,  1964]  of  Ej  (x)  ,  we  have 

ik  (z+z  )  ik  (z  +  z  ) 

co  _  o  o 

K  ,  .  e  2_  _e _ 

c  Z,Zo  ~  k  (z  +  z  )  n  k  (z+z  )(l  +  l/n) 
o  o  o  o 


iko(z  +  ZQ) 

~  rajJ:  j * - 

+  1/  k  (z  +  z  ) 


Kj(z,  z*) 


(14) 


which  agrees  with  the  approximate  solution  used  by  Jasik  [  1961]  and 
Berry  and  Chrisman  [1966].  It  corresponds  to  the  reflected  portion 
of  a  normally  incident  wave  and  is  the  first  term  of  an  asymptotic  series 
representation  intended  for  use  at  higher  frequencies  [Wait,  1962]. 

3.  DISCUSSION 

In  order  to  discuss  the  various  approximate  solutions,  the  ex¬ 
act  solution  kM,  in  (4),  which  is  expressed  only  in  terms  of  an 
integration  in  a  complex  transform  plane  a,  is  evaluated  numerically 
and  is  compared  with  our  approximate  solution  K  in  (12);  the  one 
used  by  Lavrov  and  Knyazev  [1965],  Kj^  in  (13);  and  the  one  used  by 
Jasik  [1961]  and  Berry  and  Chrisman  [1966],  Kj  in  (14).  The  earth 
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ground  is  chosen  to  have  a  conductivity  of  10  mho/m  and  a  relative 
dielectric  constant  of  10.  The  operating  frequency  is  set  at  100  MHz, 
giving  a  refractive  index  of  the  earth  n=  3. 175  +  i  0.  28  which  is  not 
particularly  large.  Both  the  amplitude  and  phase  of  these  solutions 
are  plotted  against  the  normalized  distance  k  (z  +  z  ),  i.  e. ,  the  total 
distance  of  the  source  dipole  and  the  observation  po'int,  in  Figures  1 
and  2.  When  the  source  and  the  observation  point  are  close  to  the 
ground  (i.  e. ,  kQ(z  +  z  )  s  1.0),  both  Kc  and  agree  reasonably 

well  with  the  exact  solution  id8),  while  does  not.  Perhaps 
better  than  what  we  assume  in  the  derivation,  the  error  involved  is 
less  than  10$  so  long  as  k  (zf  z}2  0,25,  or  the  total  distance  is 
longer  than  0.  04  of  a  free -space  wavelength!  Farther  away  from 
the  ground  (i.  e. ,  k  (z  +  z  )  a  1),  the  error  in  K_  is  almost  negligible, 
in  the  same  range,  however,  the  error  in  the  approximate  solution 
K  ^  tends  to  be  worse,  while  the  one  with  Kj  is  better.  This  is 
due  to  the  fact  that  the  error  in  is  on  the  order  of  l/|n2|  . 

For  all  practical  purposes,  Kc  can  be  used  adequately  for  a 
medium  with  arbitrary  refractive  index.  If,  on  the  other  hand,  |n| 
is  large  compared  with  one,  K  L  is  preferred,  because  of  its 
simplicity  for  application  to  impedance  calculations  of  low  antennas 
over  a  well  conducting  earth  [Whit,  1969]. 
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Abstract 

Hie  transmission  of  a  delta  function  plane  wave  into  a  plane  earth 
is  determined.  Hie  complex  dielectric  properties  of  the  ground 
are  modeled  using  an  empirical  relationship  that  exhibits  a  typi¬ 
cal  variation  of  conductivity  with  frequency.  Numerical  results 
are  presented  for  both  high  and  low  loss  grounds  as  a  function  of 
depth. 

Introduction 

The  subject  of  pulse  propagation  into  the  ground  has  continued  to  be  of 
general  interest  because  of  the  need  to  analyse  the  effects  on  electron¬ 
ic  systems  of  high  amplitude  electromagnetic  transients.  Such  transients 
might  be  produced  by  lightning  or  a  nuclear  burst.  The  purpose  of  this 
paper  is  to  gain  some  insight  into  the  shielding  properties  of  the  earth 
against  pulses.  Welt1  has  considered  this  problem  in  considerable  de¬ 
tail  and  gives  an  excellent  suranary  of  previous  work  by  others.  More 
recently  Kinj;  an'  Harrison^  investigated  this  problem  using  a  Gaussian- 
shaped  pulse.  The  problem  considered  here  is  that  of  a  delta  function 
plane  wave  normally  incident  on  a  plane  earth.  In  this  way,  the  re¬ 
sulting  pulse  at  different  depths  will  exhibit  features  due  to  the  lossy 
earth  which  are  not  masked  by  frequency  limitations  of  the  source  pulse. 
In  addition,  the  results  for  any  other  incident  puls-  shape  is  just  the 
convolution  of  the  incident  pulse  and  the  result  given  here.  Previous 
investigations  by  others  have  assumed  the  earth  to  have  a  conductivity 
that  is  constant  for  all  frequencies.  Hie  present  analysis  uses  a  model 
fer  conductivity  that  more  closely  matches  experimental  measurements. 

Analysis 

Hie  problem  is  solved  by  finding  the  steady  state  solution  and  then 
taking  the  inverse  transform  numerically  on  a  computer  to  find  the 
pulse  shape  at  different  depths  in  the  ground.  As  shown  in  Figure  1, 
the  '  tra^  jrmitted  wave  is  in  the  earth.  It  is  assumed  that  the 

earth  has  a  permeability  of  free  space  and  a  complex  relative  dielectric 
constant,  C3  ,  given  by 


'•o 

where  £,  ,  is  the  usual  relative  dielectric  constant  which  as  a  value 
around  8  to  25  depending  on  the  moisture  content,  a~  is  the  conductivi¬ 
ty  of  the  ground,  oj  ts  2  v  times  frequency,  and  is  the  permittivity  of 
free  space.  Normally  the  conductivity  is  assumed  constant  with  the 
frequency,  but  measurements  on  typical  soil  show  considerable  variation 
with  frequency — especially  at  higher  frequencies  where  it  tends  to  be 
directly  proportional  to  frequency,  TMs  is  well  known  for  low  loss  di¬ 
electrics  where  the  loss  tangent  defined  as  oyfcjs  e  is  usually  assumed 


independent  of  frequency.  To  account  for  this,  the  following  model  for 
I",  was  used  in  this  study: 


€j  £j 


~  J(< 


where  <rt  is  the  usual  low  frequency  conductivity  and  typically  has 
values  from  .002  to  .02  mhos  per  meter  and  L9  is  the  high  frequency 
loss  tangent.  There  is  little  information  available  on  values  of 
for  different  soils  but  measurements  on  a  soil  sample  from  the  Seattle 
area  show  a  value  of  0.5.  The  above  model  is  good  for  all  practical 
purposes  but  is  not  too  convenient  when  extended  to  the  entire  complex 
frequency  domain.  An  alternative  model  that  overcomes  this  and  still 
has  about  the  same  variation  with  real  frequencies  is 

£,  -  j  + 

The  steady  state  solution  for  the  electric  field  in  the  earth,  Et  (u>) , 
in  terms  of  the  incident  field,  £Tt  (uj)  ,  is 

E,  =  T  e; 

where  T  is  the  transmission  coefficient  of  the  earth-air  interface,  Y 
is  the  propagation  constant  of  the  earth  and  d  i3  the  depth.  T  and  Y 
are  well  known  and  can  be  written  in  terms  of  the  complex  dielectric 
constant  as 

~r  —  _ _ _ 

I  —  ,  ns? — 


jUJ  r~ 


Assuming  the  incident  pul3e  in  time  is  an  unit  amplitude  delta  function, 
then  E.  (w)  =  1.  The  field  in  the  earth  as  a  function  of  time  and  depth 


the  inverse  transform  of  E. ,  i . e . , 


£ 


,  j^t  .  /  r  e _  /  , 

e,<V'  =  rwj  Ete  =  ~J  TTJTT 

Results  -op  -*> 

The  above  transform  has  been  evaluated  numerically  on  a  computer.  The 
results  for  two  different  conductivities  and  several  different  depths 
are  shown  in  Figures  2  and  3 *  Figure  2  shows  the  results  for  0^  .ooz 

■=  10  apd  which  is  representative  of  a  moist  low  loss 

soil  found  in  the  Seattle  area.  At  all  depths,  thi  pulse  is  character¬ 
ized  by  a  positive  high  amplitude  narrow  pulse  followed  by  a  low  ampli¬ 
tude  slowly  decaying  negative  tail.  The  peak  of  the  pulse  shows  a 
delay  corresponding  to  that  expected  for  a  wave  propagating  to  the  same 
depth  in  loss  less  dielectric.  It  should  also  be  noted  that  the  posi¬ 
tive  pulse  appears  almost  symmetrical  about  its  peak.  Figure  3  shows 
the  results  for  only  a  change  in  the  low  frequency  conductivity,  i.e., 

<7^  =  .02,  =  10  and  L 9  =  .5*  This  is  representative  of  a  soil 

with  about  the  same  moisture  content  as  before,  but  with  a  greater 
amount  of  ionizing  material  present.  It  is  not  known  how  L?  might 
change  so  it  was  left  fixed  for  this  example.  The  same  pulse  character¬ 
istics  are  observed  but  with  much  lower  amplitude.  At  depths  less  than 


-t'il  _ 


2  meters,  the  positive  pulse  is  still  quite  symmetrical  and  has  a  time 
delay  the  same  as  before.  However,  at  greater  depths  the  pulse  becomes 
skewed  and  its  peak  occurs  at  a  greater  time  delay.  For  example,  at  a 
depth  of  10  meters  the  peak  occurs  at  about  250  nanoseconds  instead  of 
100  as  before.  In  addition,  the  peak  amplitude  decays  much  more  rapidly 
with  depth  than  it  did  with  the  lower  conductivity. 

Conclusions 

The  impulse  response  for  propagation  into  the  ground  has  been  found  and 
typical  examples  shown.  These  results  show  the  limiting  values  of  pulse 
width  and  rise  time  for  any  source.  For  example,  at  a  depth  of  1  meter 
the  maximum  rise  time  and  minimum  pulse  width  will  be  on  the  order  of  a 
few  nanoseconds  while  at  10  meters  depth  it  will  be  on  the  order  of  50 
to  100  nanoseconds.  The  results  for  any  other  incident  pulse  shape  can 
be  found  by  convolution  with  the  results  presented  here. 
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PERFORMANCE  OF  VALLEY  SPANNING  VLF  ANTENNAS 
BASED  ON  ELECTROMAGNETIC  MODELLING 


A.  N.  Smith  and  E.  J.  Jackson 
Westinghouse  Ge ore  search  Laboratories 
Boulder,  Colorado 


Passive  electrical  characteristics  of  numerous  transmitting 
antenna  configurations,  similar  to  valley -spanning  proto¬ 
types  of  interest  for  moderate  power  application  in  the 
OMEGA  Navigation  System,  were  measured  for  models  on 
alterable  valleys.  Resulting  parametric  curves  and 
modelled  near-zone  field  distributions  plus  soil  con¬ 
ductivities  of  candidate  sites  enable  accurate  prediction  of 
power  and  bandwidth  capability  for  cost  studies. 

The  communicator/navigator  generally  hands  his  problem  to  the 
communications  systems  engineei  in  the  form  of  a  desired  reliability 
(av.Mlability  in  some  noise  environment)  in  a  coverage  area  for  some 
signal  format  yielding  a  specified  information  rate.  These  require¬ 
ments  filter  down  to  the  antenna  system  designer  as  a  minimum 
radiated  power  (Pr)  at  some  frequency  (f)  or  range  thereof  in  some 
specified  minimum  bandwidth  (B  ).  Present-day  technology  supplies 
an  upper  limit  on  operating  voltage,  V.,  and  previous  generalized 
studies  give  a  range  of  desirable  efficiency  qag  (and  hence  input  power 
P  )  functionally  related  to  radiated  power  level  so  that  costs  will 
probably  be  near  a  minimum.  For  the  OMEGA  System,  the  respec¬ 
tive  values  have  been  established  as  P  S  10  kW,  10.  2  <  f  <  13.  6 

kHz,  B  =£  i 0  Hz  @10.2  kHz;  V.  <  25 ft  kV,  P  <  150  kW;  q  > 

,  as  i  as  'as 


At  VLF,  where  transmitting  antennas  are  almost  unavoidably  electri¬ 
cally  short,  these  six  variables  and  four  others  are  related  by  the 
four  equations  given  below  (applicable  in  general  to  any  simple  tuned 
circuit  of  lumped  parameters  and  specifically  to  the  one  illustrated): 


B 


as 


as 


,  v  .  -13  2  ,  2  2  ,4  2  _ 

=  6. 95  X  10  Cn  h  V  f  a  I  R 

0  e  t  a  r 

=  1.108  X  10"13  Cn  h2  q  _1  f4  =  f/Q 

0  e  'as  as 


=  P/i] 


as 


(MKSQ  units) 


[  1  -F(f/fQ)2]  V.  (f0  =  (ZtT  /IT  C0)_1  ) 


The  remaining  four  quantities  (electrostatic  capacity),  f  (self¬ 
resonant  frequency),  V  (mean  topload  voltage)  and  hg  (effective 
height)  are  thus  determined  if  the  original  six  are  fixed;  this  is 
usually  not  the  case,  but  certain  ones  are  allowed  to  take  on  values 
below  well-defined  upper  limits,  so  that  cost-  and  performance -trade 
comparisons  can  be  carried  out  for  several  candidate  sites  as  well  as 
for  several  configurations  in  final  single  selection.  A  significant 
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possibilities  for  tuning  helix  resistance  R^  against  ground  system 
resistance  R  ,  which  are  the  two  major  loss  components,  always  in 
regard  to  the^condition  of  antenna  system  efficiency.  In  performance 
trading,  one  finds  besides  that  specification  of  both  power  and  band¬ 
width  at  more  than  one  frequency  in  the  operating  band  over¬ 
determines  the  system.  Thus  one  specifies  minimum  acceptable  for 
each  at  that  frequency  it.  is  hardest  to  attain;  e.g.  ,  for  OMEGA 
P  at  10.  2  kHz,  and  B  (under  an  earlier  concept)  at  13. 6,  but  not 
both  at  both  frequencies  V.  must  be  corona  free. 


The  environment  is  significant  to  performance  determination  in  two 
respects:  geometrical,  relating  to  CQ,  h  (capacity  centroid),  and 
f  (distance  relative  to  wavelength  over  which  charge  musi.be  carried 
to  components  of  CQ);  and  electrical,  relating  to  the  nature  of  the 
impedance  plane,  losses  in  which  must  be  controlled  by  design  of 
some  ground  system  so  as  to  simultaneously  permit  attainment  of 
P  and  B  .  Since  the  cost  optimization  process  must  consider  both 
fixed  (inslalled)  as  well  as  operating  (capitalized  over  stated  period 
of  years),  some  latitude  is  offered  to  charge  against  site  modification 
for  efficiency  (such  as  elevating  spans  on  towers)  versus  more  trans¬ 
mitter  and  prime  power. 

The  power,  bandwidth,  voltage,  and  efficiency  figures  given  above 
imply,  through  the  fundamental  equations,  that  we  look  for  sites 
yielding  effective  heights  for  OMEGA  applications,  from  100  to  200 
meters  and  capacity  from  0.050  to  0.030  pfds;  the  inverse  relation¬ 
ship  of  capacity  and  inductance  usually  results  in  between  25  and 
40  kHz,  so  that  V  is  within  10%  of  V.  at  the  low  end  of  the  band. 
Reactances  to  tune  are  then  such  as  to  require  helix  Q's  between  1500 
and  3000  in  order  to  hold  to  values  such  that  the  residual  ground 

system  loss  R  can  be  obtained  with  reasonable  installation. 
y  g 

The  valley-spanning  approach  to  obtain  these  figures  is  advantageous 
only  because  of  the  relatively  modest  power  and  efficiency  necessary  in 
the  OMEGA  system,  since  the  very  presence  of  the  ridges  uced  as 


supports  for  the  spans  precludes  economic  installation  of  an 
extremely  elaborate  ground  system  and  also  results  in  high  field 
{both  H-  arid  E-  fields)  concentrations  in  locations  difficult  to  reach 
with  a  wire  grid.  A  high  power,  high  efficiency  installation  requires 
the  span  supports  to  be  such  as  not  to  get  in  the  way.  Displaying  the 
loss  budget  computation  for  the  ground  system  the  form 

R  =  (Rtt  +  R_) .  .  ,  +  <Rtt  +  R_)  . ,  +  R  +  R 

g  H  E  inside  H  E  outside  wires  termination 

in  which 

Rtt  i  R^  --  Re  (l/I2)  \  [ri  (HJ2  +  —  w2<  2  E2  1  dA  inside:  a.  to  a 
HE  a  J  1  0  cr  0  z  *  l  o 

outside:  a  to  1/6 
o  0 

shows  the  critical  role  that  field  distributions  play  in  making  perform¬ 
ance  predictions  have  to  do  with  loss  budget.  The  most  convenient 
way  to  arrive  at  these  is  through  scale  modelling. 

Capacity  can  be  calculated  rather  easily  to  a  fair  degree  of  accuracy, 
because  its  logarithmic  dependence  on  the  ratio  h/a  makes  this  rather 
insensitive  to  inaccuracies  in  knowledge  of  height;  not  so  with  effective 
height,  since  it  is  a  summation  of  production  of  local  height  with 
capacity.  Thus  again  modelling  by  analogue  methods  is  of  great  aid 
in  accurate  performance  predictions. 

Some  60  configurations  for  valley-spanning  antennas  for  the  OMEGA 
application  have  been  modelled  at  WGL,  at  480/1  scale  in  air.  26 
more  were  done  for  one  of  the  high-power  communication  stations. 
Another  1100  variants  were  studied  in  an  electrolytic  tank.  The 
results  generally  show  that  the  required  C.  and  h^  for  OMEGA  power 
and  efficiency  can  be  obtained  with  from  6  to  4  spans  over  valleys 
7,000  or  more  feet  wide,  with  aspect  ratio  4/1,  without  use  of 
support  towers,  and  with  ground  systems  of  the  order  of  2,000  foot 
radius  if  the  conductivity  is  10"°  mhos  per  meter  or  higher. 

Extremely  high  soil  conductivity  is  not  advantageous  if  there  is  a. 
large  seasonal  decrease  in  effective  conductivity  (e.g.  snow)  for 
which  provision  must  be  made  in  ground  system  design  to  meet 
radiated  power  requirement  because  danger  exists  that  during  high- 
conductivity  portion  of  the  year  bandwidth  cannot  be  met  without 
artifically  supplying  extra  loss.  From  field  strength  pattern  meas¬ 
urements,  no  evidence  exists  that  valley-span  antennas  of  the 
proportions  studied  exhibit  any  characteristics  of  slots,  as  was  one 
time  discussed. 

Work  performed  under  sponsorship  of  OMEGA  Project  Office,  PM-9, 
through  Chesapeake  Division,  Bureau  of  Yards  &  Docks,  U.  S.  Navy, 
under  Contract  N62477-68-C-0940. 
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SURFACE  WAVE  CONTRIBUTIONS  TO  THE  BACKSCATTERING 
FRCW  LARGE,  HIGH-DENSITY  DIELECTRIC  SPHERES. 


M.  A.  PLONUS  AND  H.  INADA 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
NORTHWESTERN  UNIVERSITY 
EVANSTON,  ILLINOIS 


An  application  of  the  Watson  transformation  to  the  Mie  series 
divides  this  solution  in  geometric  optics  and  diffracted  fields. 

It  will  be  shown  that  the  geometric  optics  fields  are  smaller  than 
heretofore  assumed.  The  diffracted  field  comes  from  two  types  of 
surface  waves.  One  of  these  is  the  dominant  contributor  to  the 
backscatter.  Numerical  results  for  values  of  ka  from  5  to  20  will 
be  presented. 


The  following  figure  shows  the  backscattering  from  a  dielectric  sphere 
as  calculated  from  the  Mie  series.  The  top  continuous  curve  shows  the  re¬ 
sults  which  are  obtained  when  the  geometrical  optics  method  is  applied  to 
the  sphere  [Atlas  and  Glover,  1962].  There  is  little  agreement  with  the 
exact  results.  After  an  application  of  the  Watson  transformation  to  the 
Mie  series  an  integral  can  be  split  off  which  can  be  identified  with  the 
geometric  optics  contribution  since  it  comes  from  the  lit  region  of  the 
sphere  [ir.ada  and  Plonus,  1969].  The  lower  curve  in  Fig.  1  show9  the  re¬ 
sults  which  were  obtained  from  this  contribution.  The  geometric  optics 
contribution  for  this  particular  range  of  ka  clearly  cannot  account  for  the 
major  backscatter. 
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Fig.  1.  Backscattering 
Cross  section  curves 
from  the  Mie  Series  and 
from  Geometric  Optics. 
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The  Watson  transformation  yields  other  terms  which  are  residue  contribu¬ 
tions  and  can  be  identified  with  surface  waves.  These  waves  can  travel  about 
the  sphere  many  times.  Fig.  2  shows  the  residue  contribution  of  one  pole 
(mode  1).  N=1  is  a  wave  that  has  travelled  around  the  back  of  the  sphere 
once,  N-2  is  a  wave  that  has  travelled  twice,  etc. 


Fig.  2.  Surface  Wave 
Contribution  to  Back- 
scattering  (dominant 
pole  only) . 


Once  the  existence  of  the  surface  waves  has  been  established,  it  is  not 
difficult  to  see  that  interference  between  waves  which  travelled  around  once, 
twice,  etc.  can  exist. 

ihe  next  figure  shows  the  interference  curves  for  all  N  for  three  separate 
poles  (1-1, 2,1).  Let  us  examine  the  resultant  curve  due  to  the  first  pole.  Up 
to  ka=7  we  have  a  simple  periodic  curve  which  when  we  re-examine  Fig.  2  could  be 
the  interference  between  N=1  and  N=2.  For  ka  >  7,  the  1=1  curve  becomes  doubly 
periodic  which  means  N=3  is  coming  into  play.  Similar  observations  can  be  made 
for  the  other  pole  contributions.  Notice  that  for  higher  l's  the  corresponding 
residue  contribution  has  a  smaller  peak  which  occurs  at  larger  ka. 


Fig.  3  Surface  Wave 
Contribution  to  Back- 
scattering  showing  inter 
ference  between  surface 
waves  of  N=l,2,3. 


The  next  figure  shows  that  the  surface  waves  are  the  major  contributors 
to  the  backscattering.  The  small  periodicity  of  ka  «  .7  can  be  explained  as 
being  due  to  interference  between  surface  waves  which  have  gone  around  the 
sphere  a  number  of  times,  whereas  the  large  periodicity  of  ka  «  10  seems  to 
be  due  to  the  contribution  of  the  different  poles  (1=1,2, ••♦)  which  peak  at 
different  values  of  ka. 


Fig.  4.  Backscattering 
Curves  showing  that  the 
surface  waves  are  the 
dominant  concributors 
to  the  total  backscatter 


H.  Inada  and  M.  A.  Plonus  (1969),  "Numerical  Results  for  Backscattering  from 
Large,  High-Density  Dielectric  Spheres",  Proc.  IEEE,  5T_,  No.  6,  1192-1193. 


D.  Atlas  and  K,  M,  Glover  (1962),  "Backscatter  by  Dielectric  Spheres  with  and 
without  Metal  Caps",  Proc.  Inter-Disciplinary  Conf.  Electromagnetic  Scattering 
(Pergamon  Press,  New  York) ,  213-236. 


-76- 


MICROWAVE  TRANSMISSION  MEASUREMENTS 
THROUGH  RADIANTLY  HEATED  ANTENNA  WINDOWS 

A.  C.  Lind,  J.  F.  Fox 

McDonnell  Douglas  Astronautics  Company  -  Eastern  Division 

and 

D.  Q,  Durant 

McDonnell  Aircraft  Company 
St.  Louis,  Missouri 


A  new  technique  has  been  developed  which  allows  continuous  microwave 
transmission  measurements  to  be  made  through  reentry  vehicle  antenna 
winders  while  they  are  subjected  to  high  heat  fluxes.  The  antenna 
window  specimen  is  radiantly  heated  by  a  cylindrical  graphite  heating 
element  which  produces  a  heat  flux  of  approximately  280  Btu/ft^-s ec . 

Unique  data  is  shown  on  the  degrading  effect  of  moisture  absorption 
in  slip-cast  fused  silica. 

1.  Introduction  -  During  reentry,  the  antenna  window  of  a  missile  is  sub¬ 
jected  to  very  high  heat  fluxes.  Communication  system  performance  is  degraded, 
not  only  by  the  plasma  that  forms  around  the  vehicle,  but  also  by  losses 
induced  in  the  window  material  as  a  result  of  the  severe  heating.  Dielectric 
constant  and  loss  tangent  data  is  generally  not  available  in  the  temperature 
region  of  interest,  that  is,  greater  than  3,000°F.  As  a  result,  others  (Mead, 
1967;  Cetaruk,  1968)  have  performed  microwave  transmission  measurements 
turough  antenna  window  samples  that  have  been  heated  with  a  plasma  jet.  This 
type  of  measurement  provides  a  qualitative  estimate  of  antenna  window  perfor¬ 
mance  during  reentry.  In  the  plasma  jet  test,  the  transmission  measurement  is 
possible  only  after  plasma  jet  shutdown  because  of  the  interference  effects 

of  the  plasma  during  the  heating  cycle.  The  interpretation  of  such  results 
can  be  misleading  if  transient  material  effects  are  important  or  if  the  window 
cools  significantly  in  the  finite  time  required  to  shut  down  the  plasma  jet. 

The  facility  described  in  this  paper  employs  radiant  heat  and  is  designed  so 
that  the  heat.'ng  element  does  not  interfere  with  the  microwave  measurement. 
Continuous  microwave  data  is  obtained  during  the  entire  heating  and  cooling 
cycles . 

2.  Description  of  Facility  -  A  schematic  of  the  facility  is  shown  in  Figure  1. 
The  antenna  window  specimen  is  supported  in  a  water  cooled,  open-ended  X-band 
waveguide  that  transmits  to  a  receiving  antenna  located  several  feet  away.  The 
waveguide  is  located  at  the  base  of  a  tubular,  electrically  heated  (45  kW), 
graphite  element  that  is  used  to  radiantly  heat  the  window  specimen.  The  sys¬ 
tem  is  purged  with  nitrogen  to  minimise  burning  of  the  graphite  element  which 
is  at  a  temperature  of  approximately  5>000°F.  A  loose  fit  between  the  specimen 
and  waveguide  results  in  poor  heat  transfer  and  minimizes  cooling  of  the  wirJow 
by  the  waveguide.  The  specimen  also  protrudes  slightly  out  of  the  waveguide  to 
insure  a  uniform  surface  temperature.  The  tubular  heating  element  concept 
allows  continuous  microwave  transmission  during  the  heating  cycle. 


rrhe  waveguide  assembly  is  interchangeable  with  an  identical  unit  containing 
a  water-cooled  copper  calorimeter  to  obtain  heat  flux  calibrations.  A  typical 
heat  flux  profile  is  shown  in  Figure  2.  A  steady  state  heat  flux  near 
280  Btu/ft^-sec  is  reached  after  approximately  30  seconds.  The  test  time  is 
arbitrary,  and  for  the  data  shown  here  was  2  minutes;  this  test  time  is  suf¬ 
ficient  to  produce  a  significant  melt  layer  on  slip-cast  fused  silica  window 
specimens. 

3.  Microwave  Transmission  Data  -  The  received  microwave  signal  is  recorded  on 
an  X-Y  recorder  which  is  calibrated  with  a  precision  attenuator  before  each 
test  in  terms  of  db  loss.  Directional  couplers  are  employed  in  the  transmis¬ 
sion  leg  to  allow  visual  monitoring  of  the  signal  generator  output  and  reflected 
power  level  during  the  test,  typical  transmission  data  at  a  frequency  of  10  GHz 
are  shown  in  Figure  3  for  slip-cast  fused  silica.  The  top  trace  was  recorded 
with  no  window  in  the  waveguide,  and  demonstrates  the  insignificant  effect  of 
the  heating  element  on  the  microwave  measurement.  Two  specimens  taken  from  the 
same  block  of  material  were  tested  at  the  heating  rate  shown  in  Figure  2.  One 
specimen  was  dried  in  a  low  temperature  oven  and  then  stored  in  a  desiccator 
prior  to  testing.  The  other  specimen  contained  a  normal  amount  of  moisture 
absorption.  The  dry  specimen  showed  a  gradual  decrease  in  signal  level  reaching 
a  steady  state  attenuation  of  0.5  db.  The  specimen  containing  the  moisture 
showed  a  sharp  decrease  in  signal  level  corresponding  to  an  attenuation  of  T  db 
and  then  gradually  returned  to  the  attenuation  level  of  the  dry  specimen. 

These  test  results  demonstrate  the  importance  of  continuously  monitoring  trans¬ 
mission  measurements  during  high  temperature  testing.  No  significant  change 
in  reflected  power  level  was  observed  during  either  test.  This  is  in  contrast 
to  the  results  of  plasma  jet  tests  where  mechanical  erosion  of  the  window  occurs, 
causing  detuning  changes  that  are  difficult  to  take  into  account  during  data 
reduction. 

U.  Summary  -  The  technique  presented  here  is  an  improvement  over  other  methods 
since  it  allows  continuous  microwave  measurements  during  the  entire  heating 
cycle.  The  data  presented  for  slip-cast  fused  silica  demonstrates  the  impor¬ 
tance  of  continuous  measurements,  and  shows  the  degrading  effect  of  moisture 
absorption  during  the  initial  heating  of  the  material,  which  could  affect 
communication  system  performance  in  a  reentry  situation.  Also,  this  technique 
does  not  produce  mechanical  erosion  of  the  window,  so  that  the  measurement  of 
window  material  losses  are  not  masked  by  detuning  effects. 
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OPTIMUM  FREQUENCY  FOR  A  VLF  TELECOMMUNICATION 
SYSTEM  USING  BURIED  ANTENNAS 


R.  Gabillard,  J.  Fontaine,  P.  Degauque 
University  of  Lille,  France 


Abstract 


This  paper  is  a  study  of  the  best  choice  of  frequency  for  a  VLF 
telecommunication  system  between  buried  receiving  and  transmitting 
antennas.  The  efficiency  of  such  a  system  depends  on  numerous 
parameters.  We  show  that  there  is  an  optimal  direction  of  the  em¬ 
bedded  antennas  which  is  different  according  to  we  use  a  frequency 
higher  or  lower  than  a  critical  one,  f  .  We  al3o  show  that  the 
electric  field  Er  increases  toward  a  maximum  value  at  a  frequency 
fGpt  •  We  give  the  expressions  of  f.c  and  fCpt  .  The  ratio  of  the 
electric  field  Ef  at  the  optimal  frequency  to  the  electric  field  E<p 
at  ELF  brings  into  evidence  the  interest  of  the  choice  of  the  optimal 
frequency  to  establish  a  VLF  subsurface  communication. 
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PROPAGATION  BETWEEN  PARALLEL  IMPEDANCE  SURFACES 
R.  B.  Dybd&l 

The  Aerospace  Corporation,  El  Segundo,  California 

The  modal  structure  of  a  parallel  plate  waveguide  characterized 
by  impedance  boundary  conditions  is  analyzed.  A  formal  solution 
for  both  parallel  and  perpendicular  modes  may  be  obtained  and 
impedance  loci  are  given  for  the  lossless  equal  impedance  case. 

The  propagation  between  parallel  impedance  surfaces  characterized  by 
impedance  boundary  conditions  is  analyzed.  This  problem  was  first 
considered  by  Barlow  and  Cullen  [1953]  and  later  by  Barlow  [1965]  and 
Wait  [1967].  The  geometry  is  depicted  in  figure  1.  The  modal  structure 


Y  JL./.J  /  f  < 
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z 


FIGURE  1.  Geometry  of  Parallel  Plate  Configuration 

may  be  formulated  from  Hertzian  potentials  [Collin,  i960]  and  a  suitable 
potential  function  is 

(1)  +H,K  -  U-P05'  -  %,s  o'Po(H  "  y)>  o-J<k2  +  P°2)iz 

The  subscript  H  goes  with  perpendicular  polarized  modes;  the  subscript 
E,  with  parallel  polarized  modes.  Real  values  of  pQ  correspond  to  bound 
surface  wave  type  modes;  imaginary  values,  to  plane  wave  components 
bouncing  between  the  surfaces. 

The  solution  for  the  parameters  p0  and  Rg,E  for  arbitrary  impedance 
values  has  been  obtained  and  the  orthogonality  relations  reveal  modal 
energies  are  coupled  in  the  lossy  case  [ Eybdal,  1968] .  For  equal  value 
impedances,  Rg  and  Rh  reduce  to  ±1  simplifying  the  modal  solution.  The 
properties  of  this  case  are  summarized  in  table  1.  Impedance  loci  may 
be  developed  for  the  lossless  case.  Note  that  losses  make  p0  complex.  The 
impedance  loci  are  depicted  in  figures  2  ana  3  and  p0  is  obtained  from 
the  intersection  of  the  loci  and  a  horizontal  line  drawn  at  the  surface 
impedance  value.  The  usual  impedance  requirement  for  bound  surface  wave 
modes  is  demonstrated  here;  i.e.  inductive  impedance  is  required  for 
hound  perpendicular  modes;  capacitive,  for  bound  parallel  modes.  The 
cutoff  condition  for  bounce  modes  is  jp0|  =  k  and  is  shown  in  the  loci 
for  a  waveguide  1.3  wavelengths  in  height. 
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TABLE  1.  Mode  Characteristics  for  the  Equal  Impedance,  Lossless  Case 


Mode  Nature  of  p0  Impedance  Value 


A 

B 

C 

D 

E 

F 

G 

H 


+1 
+  1 
-1 
-1 

+1 

+1 

-1 

-1 


Perpendicular  Modes 

real 

jMtaifcJEgiJL 

imaginary 

real 

«i£oicothl£aiS- 

W€  ? 

jkoLcot  IpoI? 

UE  2 

imaginary 

Parallel  Modes 

real 

imaginary 

real 

imaginary 

- J  cothl&J  H 

Nature  of  Impedance 


inductive 

capacitive- 

inductive 

inductive 

inductive- 

capacitive 


capacitive 

inductive- 

capacitive 

capacitive 

capacitive- 

inductive 


FIGURE  2.  Impedance  Loci  for  Perpendicular  Modes 
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IKE5DAJIC3  OF  A  DIPOL5  IN  A  LARS3  LABORATORY  rAKETOPLASMA. 

S.Y,X.  Tam 

Plasma  and  Space  Posies  Laboratory 
RCA  Limited 
Montreal,  Canada 


Impedance  of  a  3.2  GHz  ^-dipole  in  a  large  do  magnetoplasma 
(diameter  a  7  Ao  )  have  been  measured  for  various  electron  densities, 

A)  \3  A)  \2 

0$  X  =  i  — 22. \  $  2#  and  magnetic  fields,  Y»f <  1*5,  Ao  being 

V  “  /  \  /' 

the  free-spaca  wavelength,  w,  m  ,  and  w  the  signal,  plasma  and 

p0  C6 

electron  cyclotron  frequencies  respectively*  Results  have  been 
consistent  with  the  theory* 


INTEGRAL  EQUATION  APPROACH  TO  THE  RADIATION  FROM  A 
VERTICAL  MONOPOLE  OVER  AN  INHOMOGENEOUS  GROUND  PLANE 
(To  appear  in  January  1970  issue  of  Radio  Science) 

J.  R,  Wait  and  K.  P.  Spies 

ESSA  Research  Laboratories,  Boulder,  Colorado  80302 


The  possibility  that  HF  antenna  radiation  can  be  enhanced  at  low 
angles,  by  a  tapered  ground  system  is  explored.  The  model  consists 
of  a  vertical  linear  antenna  erected  over  a  radially  inhomogeneous 
ground  system.  The  surface  impedance  is  chosen  to  be  a  smooth 
monotonic  function  of  the  distance  from  the  base  of  the  antenna.  In 
a  practical  scheme,  this  can  be  achieved  by  locating  a  radial  wire 
system  in  a  symmetrical  fashion  beneath  the  monopole.  The  com¬ 
bination  of  the  equivalent  wire  grid  impedance  and  the  ground  im¬ 
pedance  is  computable  from  previously  developed  theory. 

It  is  first  shown  that  the  tangential  magnetic  field  satisfies  an 
integral  equation  of  the  Voltsrra  type.  Allowing  the  effective  surface 
impedance  to  vary  in  an  exponential  fashion,  a  power  series  solution 
for  the  tangential  magnetic  field  is  then  obtained.  ThiB  involves 
integral  and  half -integral  powers  of  the  distance  from  the  base  of  the 
antenna  to  a  point  on  the  ground  surface.  As  a  check  for  the  homo¬ 
geneous  ground  (i.  e. ,  no  wire  grid  system),  it  reduces  to  the 
appropriate  series  expansions  for  the  Somme rf eld  attenuation  function. 
The  power  series  solution,  so  obtained,  is  supplemented  by  a 
numerical  evaluation  of  the  integral  equation. 

Using  the  solutions  obtained  above  for  the  tangential  fields,  the 
radiation  pattern  of  the  antenna  is  calculated  for  several  inhomogeneous 
ground  models.  Some  earlier  results  using  simplified  assumptions  are 
confirmed.  +  It  is  shown  that  a  tapered  ground  screen  model  is  muci 
better  than  an  abruptly  truncated  system.  Some  consideration  is  also 
given  to  intentional  azimuthal  variation  of  the  ground  system. 


*J.  R.  Wait  (1?67),  Pattern  of 
ground  screen,  Can.  J.  Phys. 


a  linear  antenna  erectel  over  a  tapered 
45  (9),  3091-3101. 
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MICROWAVE  (4.75GHz)  SURFACE  IMPEDANCE  MEASUREMENTS 
OF  A  1*  DIELECTRIC  WEDGE  ON  A  PERFECT  CONDUCTOR 

R.J.  King  and  C.H.  Hustig 
University  of  Wisconsin,  Madison,  53706 


Abstract: 

When  propagating  in  the  direction  of  increasing  depth,  the 
results  near  the  wedge  apex  compare  best  with  calculations  based 
upon  simple  parallel  stratified  media  theory.  When  propagating 
in  the  direction  of  decreasing  wedge  depth,  the  results  compare 
reasonably  well  with  the  theory  given  in  1967  by  Sehlak  and  Wait. 
However,  the  nature  of  the  surface  impedance  is  decidedly  dif¬ 
ferent  for  the  two  cases,  being  inductive  near  the  wedge  apex  for 
the  first  case,  and  capacitive  with  a  negative  real  part  for  the 


second  case. 
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SESSION  VIII,  1:30  p.m. ,  Thursday,  July  17,  1969 
Round  Table  Discussion  on: 

DESIGN  TECHNIQUES  FOP.  PATTERN  CONTROL  BY  GROUND  SCREENS 
SYNOPSIS  by  R.  V.  Row,  Syl vania,  Waltham,  Mass. 


The  moderator  (R.  V.  Row)  introduced  the  two  others  sitting 
at  the  table  with  him:  Dr.  J.  R.  Wait  (ITS/ESSA  Boulder, 
Colo.),  and  R.  D.  Wengenroth  (General  Electric  HMES,  Syra¬ 
cuse,  N.Y. ),  and  proceeded  to  invite  the  participants  to 
range  further  afield  in  their  discussions  than  the  adver¬ 
tised  title  of  the  round  table  would  indicate  and  to  con¬ 
sider  terrain,  fences  and  other  obstacles  within  the  scope 
of  the  meeting. 

The  discussions  were  then  opened  by  J.  R.  Wait  who  dis¬ 
coursed  on  the  proper  use  of  transfer  impedance  and  input 
impedance  between  antennas.  He  pointed  out  that  most  radio 
propagation  work  in  the  literature  deals  only  with  the 
transfer  impedance  and  when  terms  like  transmission  loss  or 
system  loss  are  used  by  an  author  care  must  be  taken  by  the 
reader  to  ascertain  how  the  input  power  was  calculated. 
Frequently  input  power  is  calculated  on  the  basis  of  a  di¬ 
pole  in  free  space.  The  speaker  then  proceeded  to  sketch 
briefly  how  a  perturbation  technique,  often  known  as  the 
compensation  theorem,  has  been  used  to  compute  both  input 
impedance  and  transfer  impedance.  For  input  impedance  of 
an  antenna  and  ground  screen  one  need  only  be  concerned 
with  fields  close  in  to  the  antenna,  whereas  the  fields  over 
the  entire  ground  surface  lying  between  the  antennas  contri¬ 
bute  to  the  mutual  impedance.  He  commented  on  the  wide 
utility  of  the  method  and  the  fact  that  it  has  variational 
stationary  properties  which  reduce  the  relative  error  in  the 
perturbed  quantity  being  calculated,  to  a  lower  order  than 
the  error  in  the  assumed  tangential  E  or  H  fields. 

The  moderator  invited  E.  W.  Seeley  (U.S.  Naval  Weapons  Cen¬ 
ter,  Corona,  Calif.)  to  comment  on  his  calculation  of  the 
radiation  efficiency  of  a  long  horizontal  above  ground  wire 
antenna.  (See  pp.  174-178  of  Vol.  I  of  EFAP  Proceedings.) 


( 
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It  was  stated  by  Seeley  that  his  expression  for  radiation 
efficiency  is  based  on  the  actual  input  power  to  the  an¬ 
tenna.  The  moderator  attempted  to  restate  this  in  terms  of 
the  E  field  maintained  at  a  point  in  space  by  the  actual 
antenna  and  a  reference  dipole  each  with  the  same  input 
power.  Seeley  stated  that  it  is  more  useful  at  VLF  to  con¬ 
sider  power  density  radiated  in  a  certain  limited  range  of 
angles  ( 10-15  degrees)  above  the  horizon,  rather  than  total 
power  radiated,  since  this  is  the  power  that  maintains  the 
field  at  a  distance  via  ionospheric  reflection. 

The  moderator  invited  W.  L.  Curtis  (Boeing  Co.,  Seattle, 
Wash. )  to  comment  on  the  experimental  and  interpretive 
theoretical  work  done  by  he  and  Coe  [ 1964]  on  antennas  on 
large  ground  screens  over  lossy  earth.  He  pointed  out  that 
both  the  experimental  fork  and  a  subsequent  theoretical 
study  [  Curtis  1964]  using  the  compensation  theorem  showed 
that  for  engineering  purposes  the  radiation  pattern  of  any 
low  profile  antenna  on  a  finite  size  ground  plane  on  a  flat 
earth  can  be  easily  obtained  from  a  free-space  pattern  of 
the  antenna  and  ground  plane.  This  means  that  modeling  or 
other  free-space  antenna  design  techniques  can  be  used,  even 
if  the  antenna  is  to  be  placed  over  flat  ground.  To  first 
order  one  earth  loss  correction  curve  can  be  used  for  all 
frequencies  above  the  HF  band  and  all  ground  planes  greater 
than  A/2  in  radius.  The  small  change  in  this  'universal' 
correction  due  to  ground  screen  size  can  be  taken  into  ac¬ 
count  using  the  curves  in  Curtis  [ 1964] .  In  the  HF  band 
and  at  lower  frequencies  where  the  ground  refractive  index 
is  frequency  dependent  different  correction  curves  would 
need  to  be  prepared. 

Professor  L.  B.  Felsen  (Brooklyn  Polytechnic  Institute)  then 
suggested  that  there  is  another  way  of  looking  at  antennas 
and  ground  screens,  namely  ray  optics  which  provides  useful 
insight  not  only  fc  r  the  input  impedance  but  also  for  the 
radiation  field.  He  illustrated  his  remarks  by  the  example 
of  a  radiator  over  the  center  of  a  circular  screen.  In 
this  case  all  the  edge  scattered  rays  travel  in  radial  dir¬ 
ections,  and  the  scattering  properties  of  an  edge  are  well 
known  and  available  in  the  literature.  This  way  of  looking 
at  the  problem  also  provides  insight  into  the  effect  of 
non-circular  screens,  for  example  an  elliptical  screen. 
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where  only  four  rays  are  scattered  back  to  the  center. 

T.  Kaliszewski  (General  Electric  HMED,  Syracuse,  N.Y. ) 
pointed  out  that  diffraction  and  ray  optics  concepts  were 
combined  by  Senior  [ 1956]  In  studying  ground  wave  propaga¬ 
tion  across  edges.  Krause  (1967]  has  developed  similar 
Ideas  for  looking  at  sloping  perfectly  conducting  screens 
and  non-planar  screens.  Apparently  the  method  Is  not 
readily  adapted  to  Imperfectly  conducting  screens,  which  Is 
the  situation  of  greatest  practical  Interest. 

Professor  Felsen  suggested  the  the  surface  impedance  con¬ 
cept  could  be  used  with  imperfectly  conducting  screens. 

V.  Arens  (Sylvania  Electronic  Systems,  Mountain  View,  Calif.) 
commented  that  he  had  adopted  this  point  of  view  and  used  the 
plane  wave  reflection  coefficient  appropriate  to  the  specular 
point  in  question  to  calculate  the  radiation  pattern  from  a 
finite  vertical  dipole  over  an  imperfect  ground  screen. 

Edge  diffraction  was  not  considered  and  is  believed  to  be  a 
second  order  effect  [  reported  at  the  ARPA  sponsored  OHD  sym¬ 
posium,  Boulder,  Cole  1966]  .  Comparison  was  excellent 
between  predictions  m  >e  this  way  and  measurements  made  on 
a  large  circular  array  in  the  HF  band. 

J.  R.  Wait  stated  that  indeed  the  ray  approach  gives  in¬ 
sight  into  the  pattern  effects  of  ground  screens.  He  pro¬ 
ceeded  to  explain  an  approach  he  used  to  study  the  pattern 
of  an  array  of  slot  antennas  in  a  semi-infinite  ground  plane 
above  a  lossy  half  space.  Employing  the  reciprocity  principle 
he  considered  the  incident  wave  (vertically  polarized)  from 
above,  removed  the  lower  dielectric  half  space  and  replaced 
its  effect  by  an  image  plane  wave  incident  from  below.  This 
procedure  satisfies  the  boundary  conditions  except  in  a 
region  close  to  the  edge  of  the  half-screen  along  the  die¬ 
lectric  interface  and  indeed  the  results  agree  with  those 
found  by  Clemmow  [  1953]  who  used  a  somewhat  similar  plane 
wave  spectrum  approach.  The  results  also  agree  with  those 
found  from  the  compensation  theorem.  [These  results  were 
presented  at  the  Congres  International  Circuits  et  Antennes 
Hyper frequences,  Paris,  October  1957.*  Experimental  results 
dealing  with  this  same  problem  were  reported  at  this  same 

t  *  f 

J.  R.  Wait,  L*  Onde  Electrique,  38e  Annee,  No.  376  bis,  Supplement 

Special,  vol.  1,  pp.  21-29,  1957. 
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congress  by  Robieux  and  Simon.]  The  speaker  then  went  on 
to  formulate  an  exact  two  dimensional  Integral  equation 
for  Inhomogeneous  ground  problems  and  cast  in  the  form  of 
attenuation  functions.  He  also  suggested  that  the  ray 
optics  approach  is  not  easy  to  justify  in  many  practical 
cases  where  only  modest  or  even  small  size  ground  screens 
are  being  considered. 

W.  L.  Curtis  commented  that  the  ray  theory  shows  that  for 
a  truncated  ground  plane  the  radiation  pattern  has  small 
oscillatory  lobes  just  where  measurement  and  the  compensa¬ 
tion  theory  say  they  are,  and  that  this  kind  of  ’prediction’ 
works  even  on  planes  as  small  as  a  half  wavelength  in 
radius . 

The  moderator  then  invited  E.  A.  Thowless  (U.S.  Naval 
Electronics  Laboratory  Center,  San  Diego,  Calif. )  to  re¬ 
view  the  work  they  have  done  on  ground  screens.  Thowless 
briefly  reviewed  the  contents  of  four  NEL  reports  as 
follows: 

NEL  Report  1346  -  "Ground  System  Effect  on  High-frequency 
Antenna  Propagation,"  by  W.  E.  Gustavs on,  W.  M.  Chas?  and 
N.  H.  Ball!,  4  January  1966.  Based  on  model  measurements 
it  shows  that  there  is  significant  improvement  (5  to  7  dB) 
in  radiation  pattern  at  angles  below  20  degrees  due  to 
ground  screens  of  radius  2  1/2  to  7  1/2  wavelengths  and 
moderately  conducting  ground  at  frequencies  of  10  to  30  MHz. 
It  was  deduced  that  mesh  sizes  of  1  or  2  feet  would  be  ade¬ 
quate  for  screens  used  in  the  HF  band. 

NEL  Report  1359  -  "HF  Extended  Ground  Systems:  Results  of 
a  Numerical  Analysis,"  by  G.  D.  Bernard,  W.  E.  Gustavson 
and  W.  M.  Chase,  24  February  1986.  Based  on  theoretical 
work  of  Wait  and  Walters  [ 1963]  on  circular  sector  ground 
screens.  Computations  are  presented  concerning  circular 
screens  for  frequencies  of  4,  8,  16,  32  MHz,  ground  screen 
radii  of  2,  4,  8,  16,  32,  64,  128  wavelengths,  grid  spacing 
of  6,  12,  24,  48  inches  and  ground  parameters  of  e  =4, 
a  =  .003  mho/m.  €  =  10,  a  =  Q  01,  e  =  20,  a  =  .03  and  €  =  30, 
o  =  .1  and  elevation  angles  of  2,  5,  10,  15,  20,  25  degrees. 
In  general  an  improvement  of  3  to  9  dB  is  attainable  over 
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poor  soil  with  a  ground  screen  of  250  foot  radius.  The 
report  contains  a  FORTRAN  computer  program  listing. 

NEL  Report  1430,  -  "A  Numerical  Analysis  of  HF  Sector 
Ground  Screen  Systems,"  by  J.  M.  Horn,  11  January  1967. 
Presents  further  numerical  computations  based  on  Wait  and 
Walters  [  1963]  ,  this  time  for  sector  screens  (made  up  of 
radial  wires)  up  to  128  wavelengths  radius.  Frequency 
10  MHz.  e  =  10,  a  =  .01,  1A  circular  impedance  disc  of 
12  inch  square  mesh  of  No.  10  AWG  wire  for  various  sector 
angles  in  10  degree  increments  and  observing  directions 
both  on  the  symmetry  axis  of  a  sector  screen  and  at  10  de¬ 
gree  increments  away  from  this  axis.  It  was  found  that  if 
one  stays  at  least  10  degrees  away  from  the  edge  of  a  sec¬ 
tor  then  the  'gain'  at  any  elevation  angle  remains  within 
1  decibel  of  its  value  over  a  circular  ^reen  of  the  same 
radius.  It  was  pointed  out  that  the  nu^oer  of  radial  wires 
in  a  sector  was  doubled  every  time  the  screen  radius  was 
doubled.  The  use  of  such  a  radial  wire  screen  can  reduce 
the  wire  requirements  20  to  40  percent  below  a  square  mesh 
screen  having  the  same  performance.  The  report  contains  a 
FORTRAN  computer  program  listing. 

NEL  Report  1567  -  "Vertical  Plane  Patterns  of  High-frequency 
Monopoles  und  of  Elevated  Vertical  Dipoles  with  and  without 
Extended  Ground  Planes,"  by  J.  M.  Horn,  (Fall  1968).  Dis¬ 
cusses  two  methods  of  enhancing  gain  at  low  elevation  angles. 
One  considers  the  effect  of  long  ground  screens  on  a  vertical 
A/4  monopole  antenna  and  the  other  an  elevated  A/ 2  vertical 
dipole  both  with  and  without  extended  ground  planes.  Fre¬ 
quency  10  MHz  (and  a  2  to  1  frequency  range  above  and  below 
10  MHz);  ground  constants  same  as  in  NEL  report  1359.  Four 
fixed  heights  of  the  elevated  dipole  and  ground  screens  of 
4  and  8  wavelengths.  It  is  shown  that  an  elevated  vertical 
dipole  at  a  height  of  A/2  with  no  ground  plane  will  have  the 
same  low  angle  gain  as  a  A/4  monopole  over  a  4  to  8A  ground 
screen.  Apparently  more  gain  improvement  is  achieved  by 
elevating  a  dipole  than  by  constructing  a  ground  screen,  at 
least  for  the  ground  constants  examined  here.  The  impiove- 
ment  is  greater  the  greater  the  height  of  the  element.  A 
computer  program  listing  is  available  for  this  report. 
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The  moderator  invited  R.  G.  FitzGerrell  (ITS/ESSA,  Boulder, 
Colo.)  to  comment  on  a  method  for  measuring  the  absolute 
gain  of  vertically  polarized  HF  band  antennas.  FitzGerrell 
then  reviewed  his paper  [  1967]  where  he  discussed  a  method 
of  using  a  horizontally  polarized  reference  antenna  as  a 
gain  standard,  combined  with  an  airborne  platform  having 
an  antenna  whose  polarization  can  be  switched  between  ver¬ 
tical  and  horizontal.  Up  to  angles  of  about  30  degrees  the 
’gain*  of  the  horizontal  standard  is  relatively  insensitive 
(+,  1  dB)  to  the  assumed  ground  dielectric  constants  and 
conductivity  in  the  frequency  range  10  to  1000  MHz.  The 
results  of  sane  model  measurements  at  400  MHz  were  presented. 
J.  R.  Wait  asked  if  the  formula  and  graphs  presented  in¬ 
cluded  the  effect  of  ground  on  the  input  impedance  of  the 
A/2  horizontal  dipole.  The  speaker  replied  in  the  affirma¬ 
tive  and  that  his  result  agrees  with  the  calculation  by 
Surtees  f Ph.D.  thesis  at  University  of  Toronto,  1352]  . 

V.  Arens  suggested  that  on  the  airplane  shown  by  FitzGerrell 
the  close  coupling  of  the  loops  would  affect  the  actual 
radiated  power  when  one  was  rotated  from  a  nominal  vertical 
to  a  nominal  horizontal  polarization.  R.  FitzGerrell  re¬ 
plied  that  this  effect*  was  corrected  for  experimentally  in 
the  data  reduction  of  the  flight  test  measurements.  There 
was  some  further  discussion  ?mong  V.  Arens,  R.  FitzGerrell 
and  H.  Cottony  (ITS/ESSA,  Boulder,  Colo.)  concerning  im¬ 
pedance  variations  on  the  aircraft  antenna. 

The  moderator  then  drew  attention  to  the  work  done  by  A.  C. 
Wilson  [  1961]  at  CRPL  on  measuring  the  effect  of  a  ground 
screen  of  radial  wires  on  the  radiation  pattern  of  a  verti¬ 
cal  monopole. 

The  moderator  suggested  that  Professor  S.  W.  Maley  (Univer¬ 
sity  of  Colorado,  Boulder,  Colo.)  might  like  to  comment  on 
his  work  comparing  measurements  and  the  compensation  theo¬ 
rem  theory  of  ground  screens. 

Professor  Maley  referred  to  his  talk  of  the  previous  day 
in  answer  to  this  suggestion  and  invited  J.  R.  Wait  to  make 
any  additional  comments.  J.  R.  Wait  expressed  the  hope  chat 


The  actual  effect  referred  to  is  the  coupling  of  the  antennas  to  the 
aircraft  structure  (Ed.). 
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the  idea  of  a  fast  wave  screen  or  screen  which  would  exhi¬ 
bit  negative  wave  tilt  might  (mentioned  by  S.  Haley  in 
Session  IV)  lead  to  seme  interesting  new  developments. 

W.  Henry  (U.S.  Coast  Guard,  Washington,  D.  C.)  asked  how 
one  can  justify  (in  applying  the  compensation  theorem) 
keeping  the  tangential  H  fields  the  same  in  the  unperturbed 
and  perturbed  conditions.  J.  R.  Wait  stated  that  in  fact 
these  fields  aren’t  identical  but  that  due  to  the  varia¬ 
tional  (stationary)  properties  of  the  compensation  theorem 
only  second  order  errors  are  made  in  the  perturbed  quantity 
being  calculated.  If  one  wishes^  the  compensation  theorem 
can  be  used  to  verify  that  the  perturbed  tangential  H 
field  is  only  very  slightly  changed  from  its  unperturbed 
value.  The  moderator  ventured  the  opinion  that  this  sta¬ 
tionary  property  of  the  compensation  theorem  was  not  suffi¬ 
ciently  emphasired  by  its  users. 

T.  Kaliszewski  observed  that  J.  B.  Andersen  [  Teleteknik, 
Copenhagen,  1965]  measured  the  effects  of  an  extended 
ground  screen  system  on  a  log-periodic  HF  band  antenna  - 
and  obtained  excellent  agreement  with  rhe  predictions  cf 
the  compensation  theorem.*  This  agreement  should  do  much 
to  convince  those  in  doubt  about  the  credibility  of  the  com¬ 
pensation  theorem  approach.  J.  R.  Wait  by  way  of  reply 
suggested  that  it  is  really  not  necessary  to  rely  on  experi¬ 
ment  to  provide  convincing  evidence.  In  connection  with  a 
study  of  mode  conversion  of  VLF  waves  at  a  land-sea  boundary 
he  showed  [Wait,  1969]  that  a  rigorous  mode  matching  tech¬ 
nique  and  the  compens  *t ion  theorem  approach  yield  the  same 
first  order  results,  provided  the  reflected  modes  are 
neglected  in  estimating  the  fields  on  the  earth’s  surface 
in  the  later  method. 

Professor  R.  J.  King  (University  of  Wisconsin,  Madison,  Wis.) 
commented  that  the  tangential  E  field  is  oftem  more  sensitive 
to  boundary  perturbations  than  the  H  field  and  hence  esti¬ 
mates  of  the  H  field  are  likely  to  be  more  accurate  than 
those  of  the  E  field  which  is  frequently  discontinuous  at 
boundaries.  He  then  addressed  some  remarks  to  the  subject 
of  variable  surface  impedance  screens,  and  expressed  hope 
that  some  progress  will  be  made  on  the  inverse  problem  of 

J.  B ,  A  derse*.  Teleteknik  9,  No.  2,  33-40,  1965. 
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determining  from  a  specification  of  the  desired  radiation 
pattern  what  the  required  variable  surface  impedance  must 
be.  J.  R.  Wait  then  spoke  briefly  on  some  of  the  effects 
on  the  radiation  pattern  when  a  tapered  impedance  screen  is 
used.  This  way  it  is  possible  to  reduce  the  high  angle 
lobing  evident  with  abruptly  terminated  large  screens,  at 
the  expense  of  gain  at  low  angles.  The  compensation  theorem 
can  be  used  to  formulate  an  integral  equation  for  the  im- 
pedance  contrast. 

R.  J.  King  stated  that  the  solutions  to  such  a  problem  turn 
out  to  be  non-reciprocal  and  that  this  difficulty  probably 
lies  in  the  effect  of  the  local  ground  on  the  impedance  of 
the  antennas  which  must  be  somehow  properly  handled  in  the 
formulation. (JRW 8  solutions  referred  to  above  were  reciprocal. ,  Ed.) 

W.  L.  Curtis  suggested  that  the  abrupt  change  of  impedance 
at  the  edge  of  a  screen  is  necessary  to  the  attainment  of 
low  angle  gain. 

R.  D.  Wengenroth  was  then  invited  to  c raiment  on  some  of  the 
practical  problems  of  measuring  antenna  system  performance 
with  or  without  ground  screens.  He  suggested  that  even 
though  we  may  some  day  be  able  to  solve  the  synthesis  prob¬ 
lem,  it  will  be  too  difficult  to  realize  the  solution  in 
practical  terms  because  of  lack  of  control  of  the  environ¬ 
ment  over  a  great  enough  distance.  Another . difficulty  is 
the  lack  of  sufficiently  detailed  knowledge  about  the  char¬ 
acteristics  of  probing  antennas  ( in  their  environment)  used 
to  measure  performance  of  ground  based  arrays.  For  HF  com¬ 
munication  systems  needs  the  kind  of  results  currently 
available  on  measured  patterns  and  gain  are  usually  adequate. 
In  some  cases  however  better  quality  data  is  needed  and 
then  the  environmental  influences  on  the  probing  antenna 
seriously  limit  the  ability  to  make  absolute  gain  measure¬ 
ments,  but  not  relative  gain  measurements.  He  discussed 
the  need  in  many  cases  to  make  measurements  at  very  long 
distances  from  the  nominal  'antenna'  installation  to  take 
proper  account  of  hills  and  other  terrain  features  extend¬ 
ing  many  miles  away  from  a  rough  terrain  antenna  site.  In 
such  instances  a  large  high  performance  high  altitude  air¬ 
craft  like  the  KC-I35  is  needed.  The  theoretical  work  to 


-94- 


determine  what  features  of  terrain  are  influential  is  a 
necessary  concomitant  to  the  experimental  work. 

The  moderator  pointed  out  for  the  benefit  of  those  who 
did  not  attend  the  tutorial  sessions  that  R.  H.  Ott 
(ITS/ESSA,  Boulder,  Colo.)  has  developed  a  theoretically 
based  computer  executed  method  for  predicting  the  pattern 
over  irregular  terrain,  and  R.  Row  (Sylvania,  Waltham, 
Mass.)  described  a  somewhat  different  approach  to  calcula¬ 
ting  the  field  over  two  dimensional  irregular  terrain. 

V.  Arens  asked  R.  D.  Wengenroth  about  the  influence  of 
the  airborne  loop  antenna  on  cross  polarisation  measure- 
mints.  R.  D.  Wengenroth  replied  that  dome  test  results 
with  both  vertical  and  horizontal  dipoles  showed  the 
loop  to  have  ‘good1  (20-30  dB)  rejection  of  the  cross- 
polarized  component,  so  that  cross  polarization  effects 
shouldn't  have  corrupted  the  normal  measurements. 

The  moderator  thanked  the  panel  contributors  and  the 
audience  in  general  for  their  thoughtful  and  stimulating 
discussions  and  called  the  session  to  a  close  at  3:45  p.m. 
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ADDED  MATERIAL  (Ed.) 

Professor  S.  W.  Maley  suggested  during  the  round  table  discussions 
that  laboratory  models  be  used  to  study  the  effect  of  a  ’’fast-wave” 
ground  system  on  the  radiation  pattern  (with  reference  to  his  paper  on 
pg.  1-196).  However,  it  was  not  clear  how  such  a  ’’fast -wave”  im¬ 
pedance  boundary  is  to  be  fabricated. 

There  was  also  some  discussion  concerning  the  tangential  field 
approximations  used  in  applying  the  compensation  theorem  formulation 
to  the  self -impedance  calculation  of  a  monopole.  As  indicated  elsewhere, 
in  the  summarized  discussions,  the  stationary  property  permits  a  crude 
approximation  to  the  tangential  magnetic  field  to  yield  relatively  good 
self-impedance  estimates.  Professor  R.  J.  King  then  pointed  out 
formulations  involving  tangential  magnetic  fields  are  less  influenced 
by  local  inhomogeneitie-s  than  are  representations  involving  the  tangential 
electric  fields.  This  point  was  attributed  to  earlier  work  by  K.  A. 

Norton  (1942)  on  related  subjects. 
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Professor  R.  J.  King  then  suggested  that,  in  the  future,  effort 
be  directed  toward  the  synthesis  of  antenna  ground  systems.  He 
hoped  to  be  able  to  synthesize  the  system  to  meet  the  objectives, 
e.  g. ,  specify  a  variable  surface  impedance  over  the  ground  surface, 
and/or  possibly  the  ground  screen  shape  and  size. 


Reference 

Norton,  K.  A.  (1942),  The  polarization  of  downcoming  ionospheric 
radio  waves,  FCC  and  MBS  Rept.  No.  60047,  pp.  58-61  (Avail¬ 
able  through  ESSA  Library,  Boulder,  Colorado). 
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SUMMARY  OF  DISCUSSIONS  AND  POINTS  OF  VIEW  RAISED  DURING 

REGULAR  SESSIONS 


After  nearly  every  paper,  questions  from  the  floor  were  raised. 
Many  of  these  were  in  the  nature  of  comments  about  details  of  the 
presented  material  and  need  not  be  recorded  here  for  posterity. 
Furthermore,  many  of  the  points  raised,  concerning  the  effects  of 
finite  ground  plane  and  related  environmental  effects,  were  discussed 
in  the  ’’round  table”  session  on  "Pattern  Control  by  Ground  Screens.  ” 
The  discussions  and  discourses  at  this  session  have  been  admirably 
summarized  by  the  moderator  Dr.  R.  V.  Row.  His  summary  appears 
on  page  86  in  this  present  volume. 

It  was  apparent  that  many  authors  described  the  results  of  compli¬ 
cated  numerical  calculations  which  arose  from  integral  equation  or 
related  point-matching  techniques.  Since  an  analytical  view  of  the 
environmental  influence  on  antenna  radiation  leads  to  a  boundary- value 
problem,  it  is  not  really  surprising  that  various  numerical  techniques 
are  resorted  to.  This  trend  in  applied  electromagnetic  theory  has  been 
criticized  by  a  number  of  prominent  researchers  from  academia. 

Their  viewpoint  was  elegantly  put  forward  by  Professor  L.  B„  Felsen 
in  a  written  communication  following  the  c  onference.  I  am  taking  the 
liberty  to  quote  from  his  letter  as  follows: 

”If  I  may  voice  again  a  personal  concern  which  has  to  do  not 
only  with  my  session  but  with  papers  in  general,  there  is  now 
with  the  availability  of  computers  the  temptation  to  furnish  de¬ 
tailed  numerical  results  for  special  boundary  value  problems 
involving  oddly  shaped  scatterers,  propagation  through  a 
special  environment,  etc.  Such  work  is  justifiable  if  the 
numbeis  are  required  for  a  special,  particular  application. 
However,  unless  such  work  is  related  to,  and  interpreted 
quantitatively  via,  approximate  analytical  techniques  which 
grant  an  insight  into  the  radiation  and  scattering  mechanism, 
it  does  little  to  advance  our  understanding;  in  fact,  failure 
to  make  such  critical  comparisons  is  to  exploit  inadequately 
the  contribution  that  judiciously  chosen  accurate  solutions 
can  provide.  Continued  and  excessive  emphasis  on  numeri¬ 
cal  results  without  adequate  theoretical  perspective  will 
tend  to  negate  advances  made  through  the  years  in  propaga¬ 
tion  and  diffraction  theory.  I  have  often  shuddered  to  think 
of  the  day  when  orientation  toward  numerical  computation 
will  be  carried  so  far  as  to  remove  motivation  for  devices 
such  as  the  Watson  transformation,  which  not  only  do  an  in¬ 
efficient  job  much  more  efficiently  but  also  add  invaluable 
insight  into  physical  phenomena.  ” 
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The  points  raised  by  Professor  Felsen  should  be  carefully  con¬ 
sidered  by  those  involved  in  computer-oriented  approaches  to 
electromagnetics.  1  personally  share  his  misgivings  that  computer 
printout*  themselves  give  little  physical  insight.  On  the  other  hand, 
iterative  solutions  of  integral  equation  formulations  lead  to  a  quanti¬ 
tative  understanding  of  the  relevant  phenomena.  In  particular,  if  the 
variational  property  of  the  compensation  theorem  formulation  is 
utilized»an  economical  and  accurate  solution  to  a  wide  class  of  prob¬ 
lems  can  be  obtained  with  ease.  After  all,  "insight"  means  different 
things  to  different  people.  If  all  the  properties  of  a  given  electro¬ 
magnetic  environmental  problem  can  be  obtained  in  a  compact  graph¬ 
ical  chart  or  nomograph,  then  who  can  say  this  gives  less  insight 
than  an  asymptotic  result  which  has  an  uncertain  region  of  validity 
(e.  g. ,  it  may  be  only  valid  if  the  range  .r  tends  to  infinity  --  but  one 
would  like  to  know  if  ®  is  1  meter  or  )  „iegameter). 

Detailed  Comments^ 

Session  Ir.  Professor  L.  B.  Felsen  suggested  that  D.  A.  Hill  (author 
of  paper  No.  2,  pg.  1-10)  consider  the  method  of  Louis  Cagniard  in 
treating  electromagnetic  transient  problems.  (See  A.  T.  deHocp, 

Appl.  Sci.  Res.  Sec.  B,  vol.  8,  349-356,  1961.) 

Session  II:  Dr.  J.  R.  Wait  wondered  whether  R.  J.  Lytle? s  solutions 
(in  paper  No.  1,  pg.  II-  8  )  for.  the  prolate  spheroid  could  not  be  obtained 
by  applying  the  reciprocity  theorem  to  the  case  of  a  uniformly  excited 
circumferential  gap.  The  solution  of  this  scalar  problem  is  well  known. 

Several  individuals  in  the  audience  questioned  the  validity  of  Dr. 
Franc  eschetti's  results  (in  paper  No.  2,  pg.  1-33).  It  was  concluded 
that  the  results  were  only  valid  for  direct  current  or  extremely  low 
frequencies. 

C.  J.  Sletten  asked  Mrs.  Chatterjee  (paper  No.  3,  pg.  1-38)  how 
the  boundary  condition  at  the  tip  of  the  conical  structure  was  satisfied. 
Dr.  Wait  explained  that  the  eigenvalues  were  determined  from  the 
condition  that  the  tangential  electric  fields  were  zero  on  the  perfectly 
conducting  spherical  tip.  Several  questions  from  the  floor,  concerning 
the  convergence  of  the  proposed  solution,  were  left  unanswered. 


t 

The  page  number  in  Vol.  I  and  Vol.  II  where  papers  appear,  is 
indicated. 
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Several  individuals  in  t‘.ie  audience  questioned  Professor  H. 

Kurss  (paper  No.  4,  pg.  1-42)  about  the  justification  for  his  assumption 
of  uniform  current  in  the  loop.  Dr.  Wait  ventured  the  opinion  that  for 
loops  near  the  ground,  the  loop  circumference  needed  to  be  small  com¬ 
pared  with  the  skin  depth  in  the  earth.  This  's  much  more  stringent 
than  merely  requiring  that  the  loop  dimensions  be  small  compared  with 
a  free-space  wavelength.  Dr.  R.  J.  Phillips  indicated  that  lie  and  his 
colleagues  at  Berkeley  had  used  a  similar  formulation  and  the  paper 
was  published  in  a  recent  issue  of  the  journal  Geophysical  Prospecting 
(full  reference  not  available,  Ed.). 

The  interesting  paper  by  P.  R.  Bannister  (paper  No.  5,  pg.  1-45) 
provoked  considerable  discussion.  The  mathematical  validity  of  the 
image  concept  was  questioned  out  of  context  by  some  individual  in 
the  audience.  It  was  stressed  by  the  author  that  the  results  were 
valid  only  in  the  quae-i-static  sense,  i.  e. ,  where  distances  were  small 
compared  with  a  free-space  wavelength. 

C.  J.  Sletten  questioned  Dr.  D.  B.  Large  (paper  No.  6,  pg.  1-50) 
about  the  motivation  for  his  finite  wire  calculations.  The  author  re¬ 
plied  that  the  results  were  used  in  interpreting  geophysical  data  at 
ELF  for  his  company. 

Session  III:  Several  in  the  audience  asked  for  clarification  of  R.  J. 
King*s  analyses  (papers  Nos.  2  and  3,  pg.  1-64,  1-66)  concerning  wave 
propagation  over  an  impedance  boundary.  Professor  King  confirmed 
that  his  results  were  valid  even  when  the  refractive  index  of  the  ground 
was  not  large.  However,  the  transition  to  a  refractive  index  of  unity 
(i.  e. ,  free  space)  was  not  permitted.  This  appears  to  be  a  problem 
which  has  not  yet  been  solved  satisfactorily. 

Dr.  Wait  provided  a  rebuttal  (in  support  of  Arnold  Sommerfeld) 
concerning  the  severe  criticisms  of  the  "master!'  by  Carson  Tsao 
(paper  No.  6,  pg.  1-77).  In  the  rebuttal,  it  was  pointed  out  that 
Sommerfeld's  results  were  not  intended  to  be  used  in  the  quasi-static 
region  where  the  height  of  the  dipole  above  the  ground  was  small  com¬ 
pared  with  the  wavelength.  Actually,  a  quasi- static  evaluation  of 
Sommerfeld's  integrals  (published  in  1909!)  leads  to  results  fully 
consistent  with  electrostatic  concepts. 

Session  IV:  Many  of  the  points  raised  here  were  covered  in  the  round 
table  discussions.  However,  there  was  a  question  raised  by  Drs.  Row 
and  Wait  about  the  possibility  of  a  cross -polarized  radiation  from  tne 
wire  grid.  It  was  concluded  that  the  effect  is  small  when  dealing  with 
a  square  mesh  which  is  more  or  less  isotropic  in  its  reflecting 
properties.  [With  reference  to  paper  No.  6,  pg.  1-103,  by  G.  A.  Otteni.  ] 
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Se88ion  V:  In  the  papers  (Nos.  2  and  3,  pg.  1-108,  1-112)  by  K.  K. 

Mei,  et  al.  ,  the  rigidity  boundary  conditions  were  questioned.  It  was 
agreed  that  they  are  not  realistic  but  they  are  convenient. 

Several  people  asked  Dr.  Oya  (paper  No.  5,  pp.  11-39  )  how  he 
managed  to  cope  with  the  plasma  boundary  conditions  on  the  antenna 
surface.  It  appears  this  question  was  avoided  by  assuming  a  current 
distribution. 

Session  VI:  No  specific  comments  were  recorded  nor  did  the  session 
chairman  volunteer  any. 

Session  VII:  Dr.  R.  V.  Row  questioned  whether  the  anisotropic  effects 
described  by  R.  J.  Phillips  (paper  No.  1,  pg.  1-169)  would  be  present 
given  the  very  low  dc  magnetic  fields  measured  near  the  moon.  Dr. 
Phillips  assured  him  that  in  the  frequency  range  from  10-20  kHz  the 
media  would  be  anisotropic. 

E.  W.  Seeley  (paper  No.  2,  pg.  1-174)  was  criticized  for  using 
efficiency  to  characterize  the  antenna  gain  and  antenna  patterns  in 
Figures  5  and  6.  Comparison  with  established  antenna  definitions 
was  urged. 

In  G.  E.  Webber's  paper  (No.  5,  pg.  H-50),  Dr.  Wait  wondered 
about  the  appropriateness  of  calling  the  propagation  mechanism  a 
surface  wave  or  a  Zenneck  wave.  Perhaps  a  guided  wave  or  trapped 
mode  might  be  a  more  accurate  and  descriptive  term. 

Some  post  conference  interest  was  Expressed  in  the  unique  phase 
dispersion  method  described  by  P.  Cornille  (paper  No.  7,  pg.  1-190). 
Actual  data  on  a  specified  layered  geology  would  be  appreciated. 

Session  IX:  V.  R.  Arens  asked  D.  C.  Chang  (paper  No.  2,  pg.  II- 58) 
why  so  much  effort  was  made  to  obtain  closed  form  solutions  rather 
than  simply  use  a  computer  to  evaluate  the  Fourier  integral  for  all 
parameters  of  interest.  This  was  answered  by  Dr.  Chang  with  the 
comment  that  a  closed  form  solution  is  preferable  because  of  the 
greater  insight  into  the  solution  which  it  provides  and  because  of  its 
greater  utility  in  design.  Dr.  Wait  (co-author)  also  commented  later 
that  the  closed  form  solutions  were  needed  in  the  nvtmerical  inversion 
of  the  integral  equations  for  a  finite  linear  antenna  over  a  conducting 
half-space.  Also,  there  is  an  economic  aspect  of  the  problem  which 
is  not  negligible  unless,  of  course,  one  has  unlimited  funds. 

Comments  on  the  physical  realizability  of  the  model  used  by  Dr. 
Curtis  (paper  No.  3,  pg.  11-64  )  were  made  by  Dr.  Wait.  He  indicated 
that  the  postulated  complex  dielectric  constant  did  not  satisfy  the 
Kramer s-Kronig  relations  and  thus  apparent  violations  of  causality 
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could  be  expected.  Dr.  Wait  and  Dr.  Curtis  agreed,  however,  that 
the  violation  was  not  practically  very  important. 

A  comment  was  made  by  Pitt  Arnold  after  the  session  that,  in 
his  experimental  investigations,  he  has  never  observed  minima  in 
the  vertical  radiation  patterns  like  those  shown  by  Dr.  J.  R.  Wait 
(paper  No.  10,  pg.  11-84  ).  Mr.  Arnold  further  commented  that  in 
experimenting  with  various  types  of  terminations  for  radial  wire  ground 
systems  the  radiation  pattern  was  negligibly  affected.  Professor  S.  W. 
Maley  commented  that  differences  were  probably  due  to  different  ground 
and  ground  system  parameters.  Dr.  Wait  later  commen'ed  that  the 
calculated  lobe  structures  for  truncated  ground  systems  have  been 
verified  by  W.  E.  Gustafson  and  colleagues  at  the  Navy  Electronics 
Laboratory  and  W.  L.  Curtis  and  colleagues  at  Boeing  Laboratories, 
Seattle.  Of  course,  for  a  tapered  system  (including  some  radial  wire 
systems),  the  lobe  structure  is  less  pronounced.  Dr.  D.  B.  Large 
asked  how  an  approximate  solution  having  no  mathematical  basis  can 
be  verified  as  meaningful.  Dr.  Wait  answered  that  certain  aspects 
of  a  solution  such  as  boundary  behavior  or  asymptotic  behavior  can  be 
examined  for  correctness  or  approximate  correctness.  Also  mathe¬ 
matical  means  are  available  for  verifying  that  an  approximate  solution 
is  nearly  equal  to  the  exact  solution  of  the  equation  in  question. 

A  question  concerning  R.  J.  King’s  paper  (No.  11,  pg.  11-85  ) 
was  raised  about  the  measurement  technique  for  the  electric  field 
over  the  dielectric  wedge.  Professor  King  explained  that  measure¬ 
ments  were  made  by  a  scattering  dipole  technique  and  that  great 
difficulty  was  experienced  in  synthesizing  a  low  reflection  termination 
at  the  edges  of  the  dielectric  wedge  material. 

Mrs.  G.  B.  Goe  of  ITS/ESSA  made  an  unscheduled  presentation 
on  Noise  and  Cerenkov  Radiation.  (See  pg.  11-114)  ' 


1  have  prepared  this  summary  of  the  discussions  with  the  help  of 
written  comments  supplied  to  me  from  C.  J.  Sletten  and  S.  W.  Maley 
for  sessions  which  they  -haired.  All  the  other  comments,  except  for 
the  quoted  gem  from  Professor  Felsen  are  based  on  my  own  (incomplete) 
recollections  of  the  events.  I  accept  the  blame  for  any  baises  that  the 
reader  may  detect.  Certainly,  they  do  not  reflect  the  official  (or  un¬ 
official)  policy  of  the  federal  agencies  of  the  U.  S.  Government  which 
have  sponsored  this  meeting. 


James  R.  Wait 
4  August  1969 
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CORRECTIONS  TO  VOLUME  I 


Finite  tubular  antenna  above  a  conducting  half  space  (Paper  No.  1  in 
Session  HI.  pp.  59-63)  by  D.  C.  Chang 


Because  of  errors  involved  in  numerical  computation,  the  scale 
factor  on  the  abscissae  in  Figs.  2  and  3  should  be  modified  by  about  a 
factor  of  2.  However,  the  discussion  following  these  figures  is  still 
valid.  The  correct  input  conductance  (in  milli-mhos)  of  a  tubular  an¬ 


tenna  of  radius  a  =  0.  007  X.  and 


h  =  0.  25  is  tabulated  in  the  following: 


Ci 


0.27 


0.3 


0.  5 


1.0 


2.0 


4.0 


Sea -water 
(o  =4  ohm/m, 

er  =  SO  at  600  MHz)  7.  14  8.  32  9. 09  8.  88  8.  84  8.  84 

Wet-earth 

(a  =  10  X  10“ 3  mho/m, 

er  =  10  at  100  MHz)  7.66  8.63  8.97  8.  86  8.  84  8.84 


When  the  antenna  is  located  several  wavelengths  above  the  ground, 
the  theoretical  results  compare  favorably  with  experimental  measure¬ 
ment  of  a  single  antenna  in  the  free-space^  (8.90  miPi-mhos)^  For  a 
more  detailed  discussion,  the  reader  is  referred  to  a  paper  by  Chang 
and  Wait  (1969)  listed  in  the  references. 

t  R.  B.  Mack,  "A  study  of  Circular  Arrays,  ”  Tech.  Rept.  No.  382, 
Crufts  Lab. ,  Harvard  Univ. ,  Cambridge,  Mass.,  May  1,  1963. 


Reflection  of  waves  of  arbitrary  polarization  from  a  rectangular  mesh 
ground  screen  (Paper  No.  6  in  Session  IV,  pp.  103-107) by  G.  A.  Otteni 

On  page  104,  the  last  paragraph  should  read:  It  is  assumed  that. ...  .as 
E  =  Z.  I,  (Z.  is  the  internal  impedance  of  the  wire  and  is  approximately 
Z  /Zit  r)  .  „ . .  for  an  axial  current  I  . . .  .with  surface  impedance  Z  . 
In'addition. . . 

On  page  107,  for  the  sample  calculations  in  the  second  paragraph,  the 
grid  wire  spacing  |a  =  »  and  b  =  2  m.  for  parallel  wires. 


Numerical  solution  of  dipole  radiation  in  a  compressible  plasma  with  a 
vacuum  sheath  surrounding  the  antenna  (Paper  No.  3  in  Session  V, pp. 
112  -116  by  S.  H.  Lin  and  K.  K,  Mei 

The  second  sent  nee  of  the  last  paragraph  on  page  114  should  have 
read:  This  behavior  of  input  reactance  near  0?  is  in  contrast  with  that 
for  infinitely  long  cylindrical  antenna  found  by^Miller  (1968). 
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Phase  measurements  of  electromagnetic  field  components  (Paper  No. 
7  in  Session  VH,  pp.  190-195)  by  P.  Cornille 

The  corrected  version  of  Fig.  4  is  given  below: 


-104- 


Radiation  from  a  parallel-pjate  waveguide  into  an  inhornogeneously 
filled  space  {Paper  No.  6  in  Session  t,  pp.  28-3})  by  R.  J. 

Kostelnicek  and  R.  Mittra  -  revised  pg.  29 

for  the  canonical  problem.  The  first  step  involves  expressing  the 
fields  in  region  B  in  terms  of  waveguide  modes  and  in  region  E  in  terms 
of  a  continuous  spectrum  (Fourier  integral)  representation.  Similar 
transform  representations  are  employed  in  regions  A,  D,  and  C.  Typical 
representations  for  the  magnetic  intensity  in  regions  B  and  E  for  a 
TM^  incident  mode  of  unit  amplitude  are  given  by 

00 

(Hy)B  -  cosj^Je  +  ]T  Bg  cos(^)  e  8  (1) 

s«o 


-  / 


-5(z-£-2t) 

E(a)e  cos (ax) da 


The  longitudinal  wave  numbers  are  given  by  6  =  V (st/b) 2-k2  in 

_ s  o 


/  2  T  /  2  b 

sgion  B,  £  -  a  -k  in  E,  and  n  *»  VY  -k*' 


k  in  E,  and  n  =  VY  -k~  in  C.  The  mode  coefficients 
o  o 


and  spectral  weight  functions  are  B  ,  D(a) ,  E(a),  C(y),  and  A (a). 

s 


The  next  step  is  to  solve  for  the  various  mode  and  weight  co¬ 
efficients  by  matching  the  transverse  field  components  at  z  =  £  and 
z  =  H  +  2t,  thereby  obtaining  D(a)  and  A (a)  in  terms  of  E(a).  Field 
matching  is  once  more  carried  out  at  the  z  =  0  interface  for  |x| 
and  subsequently  for  |x|  ^  b.  The  resulting  four  equations  are  Fourier 
transformed  and  combined  to  yield  the  following  relationships: 


b(l+6°)5° 
P  s 


(3) 


00 

cos  (Yb)TfT)E (Y)R^(y)  =  P.  V.  J  E'(a)  f  -j;  -  f--- da,  (4) 

0  c  '  n  5  +  nj 

together  with  two  companion  relationships  giving  the  coefficients  B 
and  weight  function  C(Y)  in  terms  of  an  integral  relationship  in-  s 
volving  E(a).  In  equations  (2)  and  (3)  we  have  written 
E'(a)  =  a  sin(ab)  E(a),  R'(a)  =  R(a)  exp  (£2,)  and  Q"(a)  =  Q(a)  exp(-£Jl). 
The  quotients  1/R(a)  and  Q(a)/R(a)are  respectively  the  transmission 
and  reflection  coefficients  for  a  uniform  plane  wave  incident  onto 
the  dielectric  cr  plasma  slab  at  an  angle  given  by  9  =  sin  (a/k  ) . 
Equations  (3)  and  (4)  represent  a  homogeneous  integral  equation  ° 
for  E(a)  and  a  requirement  that  E(a)  take  on  certain  specified  values 
at  a  =  sTT/b  (s  =  0,  1,  2,  . ). 
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The  above  f< — collation  is  exact,  and  for  the  limiting  case  when 
K  =  1,  an  exact  solution  for  E(a)  can  be  constructed  using  the 
function-theoretic  technique. 

Modified  Method  of  Solution 

The  solution  of  (3)  and  (4)  for  E(ct)  by  the  modified  function- 
theoretic  technique  is  accomplished  by  the  construction  of  a  sectionally 
holomorphic  function  F(<u)  of  a  complex  variable  to,  which  has  a  certain 
pole-zero  configuration,  and  specified  branch  singularities. 

Integrations  in  the  complex  plane  yield  results  which,  when  compared 
with  (3)  and  (4)  and  the  compaion  relations  for  B  and  C(y),  give 
the  required  solutions  and  the  normalization  condition.  These  are: 


F(-8  )  =  -b(-l)S8  (1+6°)B  (5) 

s  s  s  s 

F(5)  =  ir^R"(a)ejabE(a)  (6) 

f(i)r^(y1  +  F(-n)  =  'rmC(Y)  (7) 

and 

F (8  )  =  b(-l)P8  (1+6°)  (8) 

P  P  P 


The  function  F(u>)  may  be  factored  into  the  form  F(w)  =  F^Uo)  T(io)  , 
where  F^  (to)  represents  that  function  employed  in  the  solution  of  the 
canonical  problem  (Mittra  and  Bates,  1965)  ,  and  T(U))  represents  the 
departure  from  the  canonical  function  due  to  the  presence  of  the  slab. 
The  preceeding  integrations  and  comparisons  together  with  the  factored 
form  of  F(u)  yield  an  auxiliary  integral  relationship 

T "  1  +  f  A(z)  T (z)  dz  .  (9) 

u-8  w  +  z 

P  a 

Equation  (9)  is  not  very  convenient  for  numerical  methods  of  solution 
when  a  is  on  the  path  0  due  to  the  singular  nature  of  the  partial 
kemtl  A  (z) ,  which  has  poles  on  a  due  to  the  surface  modes  excited 
with:  i  the  dielectric  slab.  However,  the  path  a  may  be  deformed  to 
say  o'  on  which  the  integrand  of  (9)  is  wholly  analytic.  Numerical 
methods  are  now  employed  together  with  a  process  of  analytic  continu¬ 
ation,  and  the  required  values  T(w)  and  T (-to)  obtained. 

The  near  fields  may  be  obtained  from  the  modal  expansion  in 
region  B  and  from  the  Fourier  transforms  in  the  open  regions.  The 
radiation  fields  are  obtained  directly  from  the  spectral  weight 
coefficients  by  employing  the  method  of  saddle-point  integration 
(Collin,  1960).  The  surface  modes  nay  be  obtained  from  the  residues 
of  the  integrand  in  the  transform  representation  in  region  A. 
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NOISE  AND  CERENKOV  R >  DIATION 


G.  B.  Goe 

Institute  for  Telecommunication  Sciences 
Environmental  Science  Services  Administration 
Boulder,  Colorado  80302 

Abstract 

Radiation  generated  by  a  moving  electron  is 
suggested  as  a  source  of  noise  in  the  ionosphere. 

The  moving  electron  travels  with  a  velocity  which 
exceeds  both  the  thermal  velocity  of  the  electrons 
in  the  ambient  (homogeneous)  plasma  and  the  phase 
velocity  of  light  in  the  medium  in  order  to  fulfill  the 
Cerenkov  condition. 

Introduction 

It  has  been  convincingly  demonstrated  by  using  automatic  gain 
control  (AGC)  records  that  the  principal  plasma  resonances  observed 
by  the  Alouette  I  topside  sounder  do  not  occur  naturally  but  are  instead 
excited  by  the  sounder  transmitter.  A  "noise"  background  is  recorded, 
however,  even  when  the  transmitter  is  turned  off.  Theoretical  re¬ 
sults  suggest  that  the  Cerenkov  mechanism  is  useful  to  interpret  the 
presence  of  the  background  noise  and  subsequent  noise  cut-off  when 
the  transmitter  is  turned  off  or  the  appearance  of  a  resonance  spike 
when  the  transmitter  is  turned  on  and  the  antenna  is  providing  energy 
to  the  electrons. 

Theory 

Dispersion  relations  are  derived  using  kinetic  theory.  The 
Cerenkov  coherency  condition  is  formally  introduced  into  newly  de¬ 
rived  plasma  dispersion  relations  in  order  to  arrive  at  algebraic 
expressions  for  the  dispersion  relationships  of  a  fully  ionized 
warm,  magneto-active  plasma  in  the  presence  of  a  moving  charge. 
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Comments 

The  particular  restriction  that  the  moving  electron  be  con¬ 
strained  to  travel  parallel  to  the  earth's  magnetic  field  gives  rise 
to  a  well-defined  cut-o^f  for  the  Cerenkov  zone  at  the  upper  hybrid 
frequency  (f^  +  f|j  )a  where  f  jyj  is  the  plasma  frequency  and 
f  is  the  cyclotron  frequency  of  the  electrons.  Under  the  same 
circumstances,  the  conditions  near  the  plasma  frequency  are  ill 
defined  and,  among  other  things,  temperature  dependent. 

This  theoretical  work  was  published  in  detail  in  a  research 
report:  Goe,  G.  B.  ,  "Plasma  Resonance  and  Radiation  in  the  Upper 
Ionosphere,  ”  School  of  Pharmacy,  University  of  Colorado,  28  July 
1967,  and  has  been  out  of  print  for  some  time.  Additional  copies  of 
the  report  are,  however,  now  available  from  the  Department  of 
Electrical  Engineering,  University  of  Colorado,  Boulder,  Colorado. 
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